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ABSTRACT 

We present new near-infrared J and K imaging data for 67 galaxies from the 
Universidad Complutense de Madrid survey used to determine the SFR density of the 
local universe by Gallego et al. (1995). This is a sample of local star-forming galaxies 
with redshift lower than 0.045, and they constitute a representative subsample of the 
galaxies in the complete UCM survey. From the new data, complemented with our 
own Gunn-r images and long-slit optical spectroscopy, we have measured integrated 
-fT-band luminosities, r — J and J — K colours, and Ha luminosities and equivalent 
widths. Using a maximum likelihood estimator and a complete set of evolutionary 
synthesis models, these observations have allowed us to estimate the strength of the 
current (or most recent) burst of star formation, its age, the star-formation rate and 
the total stellar mass of the galaxies. An average galaxy in the sample has a stel- 
lar mass of 5x10 10 Mq and is undergoing (or recently completed) a burst of star 
formation involving about 2 per cent of its total stellar mass. We have identified two 
separate classes of star- forming galaxies in the UCM sample: low luminosity, high exci- 
tation galaxies (Hll-like) and relatively luminous spirals galaxies (star burst disk- like). 
The former show higher specific star formation rates (SFR per unit mass) and burst 
strengths, and lower stellar masses than the latter. With regard to their specific star 
formation rates, the UCM galaxies are intermediate objects between normal quiescent 
spirals and the most extreme HII galaxies. 
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1 INTRODUCTION 

The study of the evolution of the Star Formation Rate (SFR) 
of individual galaxies and the SFR history of the Universe 
has experienced considerable progress recently (see, e.g., 
Madau, Dickinson & Pozzeti 1998 and references therein). 
These are key observables needed to extend our understand- 
ing of galaxy formation and evolution. In the last few years, 
the combination of very deep ground-based and HST multi- 
band imaging with deep spectroscopic surveys carried out 
with 4-m and 10-m class telescopes has allowed the sketch- 
ing of the SFR history of the Universe up to z > 4 (see, e.g., 
Madau et al. 1998 and references therein). 

A great deal of effort has been devoted to both obser- 
vational and theoretical studies of star-forming objects and 
their evolution with look-back-time. Deep imaging and spec- 
troscopy of faint galaxies at intermediate and high redhsifts 



have yielded vast amounts of quantitative information in this 
field (Lilly et al. 1995; 1998 and references therein; Driver, 
Windhorst & Griffiths 1995; Steidel et al. 1996; Lowenthal et 
al. 1997; Hammer et al. 1997; Hu, Cowie & McMahon 1998; 
see Ellis 1997 for a recent comprehensive review). Although 
substantial uncertainties still exist, a reasonably coherent 
picture is emerging. The Star Formation Rate density of the 
universe was probably about an order of magnitude higher 
in the past than it is now, perhaps peaking at z ~ 1-2 (e.g., 
Gallego et al. 1995; Madau et al. 1996; Connolly et al. 1997; 
Madau et al. 1998). These observational results seem to be 
in good agreement with the predictions of recent theoreti- 
cal models of galaxy formation (Pei & Fall 1995; Baugh et 
al. 1998; Somerville, Primack & Faber 1999), although the 
question of whether the SFR density decreased beyond z ~ 2 
is still a matter of intense debate (Hu et al. 1998; Somerville 
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et al. 1999; Hughes et al. 1998; Barger et al. 1998; Steidel et 
al. 1999). 

Given the large redshift range covered by these stud- 
ies, different SFR indicators have perforce been used, all of 
which have different calibrations, selection effects and sys- 
tematic uncertainties. These indicators include emission line 
luminosities (e.g., Ha, H/3, [OIIJA3727A), blue and ultravio- 
let fluxes, far-infrared and sub- mm fluxes, etc (see, e.g., Gal- 
lego et al. 1995; Rowan- Robinson et al. 1997; Tresse & Mad- 
dox 1998; Glazebrook et al. 1999; Treyer et al. 1998; Madau 
et al. 1996; Connolly et al. 1997; Hughes et al. 1998; Barger 
et al. 1998; see also Chariot 1998 and Kennicutt 1992). It 
would be highly desirable to use the same SFR indicator at 
all redshifts, so that the problems related to different selec- 
tion effects and systematics could be avoided. It is widely 
accepted that the Ha is one of the most reliable measure- 
ments of the current star formation rate (modulo the IMF; 
see, e.g., Kennicutt 1992). Several groups have used the Ha 
line to estimate SFRs at different redshifts, from the local 
universe to beyond z = 1 (Gallego et al. 1995; Tresse & 
Maddox 1998; Glazebrook et al. 1999), albeit with very dif- 
ferent sample selection methods. Nevertheless, it is clear that 
it is now necesary to build sizeable samples of Ha-selected 
star-forming galaxies at different redshifts and study their 
properties. One would like to know the preferred sites of star 
formation in the local universe and beyond, and the main 
propeties of the star-forming galaxies and their evolution. 
Some questions that need to be answered include: does star 
formation mainly occur in dwarf, starbursting galaxies or in 
more quiescent, normal L* galaxies? how has that evolved 
with time? what fraction of the stellar mass of the galaxies 
is being built by their current star-formation episodes? 

Progress towards answering questions such as these re- 
quires, as a first step, a comprehensive study of the proper- 
ties of the star-forming galaxies in the local universe. The 
Universidad Complutense de Madrid survey (UCM here- 
after; Zamorano et al. 1994, 1996) is currently the most com- 
plete local sample of galaxies selected by their Ha emission 
(see section ^). It has been used to determine the local Ha 
luminosity function, the SFR function and the SFR density 
(Gallego et al. 1995). It is also widely used as a benchmark 
for high redshift studies (e.g., Madau et al. 1998 and refer- 
ences therein). Thus, the UCM survey provides a suitable 
sample of local star-forming galaxies for detailed studies. 

Both optical imaging (Gunn-r; Vitores et al. 1996a, 
1996b) and spectroscopy of the whole UCM sample (Gallego 
et al. 1996; GAL96 hereafter; see also Gallego et al. 1997) 
are already available. The optical data provides information 
on the current star-formation activity, but is rather insensi- 
tive to the past star-formation history of the galaxies. In this 
paper we present new near-infrared imaging observations for 
a representative subsample of UCM galaxies. The near in- 
frared luminosities are sensitive to the mass in older stars, 
and therefore provide a measurement of the integrated past 
star formation in the galaxies and their total stellar masses 
(see, e.g., Aragon-Salamanca et al. 1993; Alonso-Herrero et 
al. 1996; Chariot 1998). Alonso-Herrero et al. (1996; AH96 
hereafter) carried out a pilot study of similar nature with 
a very small sample. We will now extend the work to a 
galaxy sample that is large enough for statistical studies, 
and that is expected to represent the properties of the com- 
plete UCM sample and thus those of the local star-forming 



galaxy population. We will also improve the work of AH96 in 
two fronts: first, we will use up-to-date population synthesis 
models, and second, we will use a more sofisticated statisti- 
cal technique when comparing observational data and model 
predictions. 

In section ^ we briefly introduce the UCM sample. In 
section |^ the observations, reduction procedures, and data 
analysis are described. The evolutionary synthesis models 
are presented in section 0, and the results are described in 
section ^[ Finally, section g contains a summary of this work. 



2 UCM SURVEY 

The UCM survey is a wide-field objective-prism search for 
star-forming galaxies which used the Ha emission line as 
main selection criterium (Zamorano et al. 1994, 1996). This 
survey was carried out at the 80-120cm Schmidt Telescope 
of the Calar Alto German- Spanish Observatory (Almeria, 
Spain), using IIIaF photographic plates. The identification 
of the emission-line objects was done by visual inspection 
of the plates over the 471.4 square degrees that the survey 
covers. An automatic procedure for the detection has also 
been developed by Alonso et al. (1995, 1999) which avoids 
possible human subjectivities in the selection. The number 
of emission-line candidates found was 264, about 44 per cent 
of them previously uncatalogued. This yielde d a detection 
rate of about 0.6 objects per square degree (Zamorano et 



al. 1994). A total of 191 of these objects were confirmed 



spectroscopically by GAL96 as emission-line galaxies. 

The wavelength cut-off of the photographic emulsion 
limits the redshift range spanned by the survey to Ha- 
emitting objects below z=0.045±0.005. The completeness 
tests performed (Vitores 1994; Gallego 1995) ensure that 
the Gunn-r limiting magnitude of the whole sample is 16. 5 m 
with an Ha equivalent width detection limit of about 20A. 

Details about the observations, data reduction, reliabil- 
ity and accessibility of the complete data set are summarized 
in Zamorano et al. (1994, 1996). 



3 OBSERVATIONAL DATA 



3.1 Optical imaging 



The comple te description of the optical Gunn-r ( Thuan & 
Gunn 1976) observations and image reduction is given in 
Vitores et al. (1996a, 1996b). Briefly, these images were ac- 
quired during a total of eight observing runs from December 
1988 through January 1992 using different CCD detectors on 
the CAHA/MPIA 2.2-m and 3.5-m telescopes, both at Calar 
Alto (Almeria, Spain). 

3.2 Near-infrared images 

Near-infrared (nIR hereafter) images in the J (1.2^im) and 
K (2.2/xm) or Tf'-bands (2.1/xm), were obtained for 67 galax- 
ies from the UCM survey during three observing runs at the 
Lick Observatory and one run at the Calar Alto Observa- 
tory. 

The three Lick observing runs took place in 1996 
(January 9-14, May 4-7 and June 7-9). We used the 
Lick InfraRed Camera (LIRC-II) equipped with a NICMOS3 
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256x256 detector on the 1-m telescope of the Lick Observa- 
tory (California, USA). The instrumental setup provided a 
total field of view of 2.4x2.4 square arc minutes with a spa- 
tial scale of 0.57" per pixel. For details about the LIRC-II 
camera and the NICMOS3 detector used see Misch, Gilmore & 



Rank (1995). The Calar Alto observations were carried out 



in 1996 (August 4-6). We used the MAGIC camera with a 
NICMOS3 256x256 detector attached to the CAHA/MPIA 
2.2-m telescope at Calar Alto (Almeria, Spain). The field of 
view was 2.70x2.70 square arc minutes and the spatial scale 
0.63" per pixel. Details abou t the MAGIC camera can be 
found in Herbst et al. (1993). Images were obtained in the 
J and K' bands in all the observing runs, except for the 
Janua ry 9-12 one, when a stand ard K filter was used. A K' 
filter (Wainscoat & Cowie 1992) was used in order to reduce 
the thermal background introduced by the red wing of the 
standard K passband. 

The observational procedure fol lowed was extensively 
described in Aragon-Salamanca et al. (1993). Briefly, we sub- 
divided the total exposure time required for each object in 
a number of images, offset by a few arcseconds, in order 
to avoid saturation. Also, blank sky images were obtained 
between consecutive object images for sky-subtraction and 
flat-fielding purposes with offests larger than ~ 1 arc minute. 
This procedure allows to reduce the effect of pixel-to-pixel 
variations, bad pixels, cosmic rays, and faint star images in 
the sky frames. 

The reduction was carried out using our own IRAlf] 
procedures, following standard redu ction s steps also de- 
scribed in Aragon-Salamanca et al. (1993), these included 
bad pixel removal, dark substraction, flat-fielding, and sky 
substraction. Finally, all the object images were aligned, 
combined and flux calibrated. Flux calibration was per- 



formed usi ng sta ndard stars from the lists of Elias ( 1982 ) and 
Courteau (1995) observed at airmasses close to those of the 



objects. The atmospheric extinction coefficients used were 
kj=0.102 mag/airmass and kk =0.09 mag/airmass, while in- 
dependent zero-points were derived for each night. The K'- 
band magnitudes were converted into K-band magnitudes 
using the empirical relation given by Wainscoat & Cowie 
( [l992| ), K' - K=0.22x(H - K). B ased on the zero-redshift 
SEDs of Aragon-Salamanca et al. ( 1993] ) and the mean nIR 
colours given in AH96 for a small sample of UCM galaxies, 
we used an H - K colour of 0.3±0.1 m . Thus, the A" - K 
correction was 0.07±0.02 m . 

Aperture photometry was carried out on the Gunn-r 
and nIR images using the IRAF/APPHOT routines. We 
measured rJK magnitudes and optical-nIR colours (r — J, 
J — K) through several physical apertures, including three 
disk-scale lengths^ (see Vitores et al. 1996a). We also mea- 
sured total Tf-band magnitudes using physical apertures 
large enough to ensure that all the light from the galaxy was 
included. Finally, we corrected for contamination from field 
stars by replacing affected pixels by adjacent sky counts. Cir- 



* IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which is operated by the Association of Universities for 
Research in Astronomy, Inc. (AURA) under cooperative agree- 
ment with the National Science Foundation. 

t _ffo=50kms -1 Mpc - 1 and <jo=0.5 have been assumed through- 
out this paper 



cular aperture r — J and J — K colours measured at three 
disk-scale lengths and integrated 7f-band magnitudes are 
given in Table [lj Magnitude and colour errors include both 
calibration and photometric uncertainties. In all the ob- 
jects analyzed, except UCM0014+1748 and UCM1432+2645 
(which was observed off-centre), the field covered by the de- 
tector was large enough to include three disk-scale length 
apertures. For these two galaxies we obtained the three disk- 
scale colours from their extrapolated growth curves. The dif- 
ferences between the larger measurable aperture and the ex- 
trapolated values were 0.05 m for UCM0014+1748 and 0.1 m 
for UCM1432+2645. 

In our analysis, the rJK magnitudes and optical- 
nlR colours are corrected for Galactic and internal ex- 
tinction. Since most of the luminosity of these galaxies 
in the rJK passbands comes from the stellar continuum, 
we have estimated the colour excesses on the continuum, 
E(B — V) continuum, with the expression given by Calzetti, 
Kinney & Storchi-Bergmann (1996, see also Calzetti 1997a; 
Storchi-Bergmann, Calzetti & Kinney 1994), 



E(B - V) 



continuum 



= 0.44 x E(B - V) s 



(1) 



where the E(B — V) sas were obtained from the spectroscopic 
Balmer decrements measured by G AL96 (see below) . We as- 
sumed a diffuse dust model that implies a total-to-selective 
extinction ratio of i?y-=3.1 (see Mathis 1990; Cardelli, Clay- 
ton & Mathis 1989). Assuming a Galactic extinction curve, 
we obtained that A r /A y, Aj/Ay an d Ak/Av are 0.83, 0.28 
and 0.11 respectively (Vlathis 1990). Note than when cor- 
recting the Ha fluxes and equivalent widths we assume that 
the line emission comes from the gaseous component, but 
the continuum is mainly stellar. In table [l] we present the 
observational data before correcting for extinction, together 
with the gas colour excesses needed for the correction. 

3.3 Optical spectroscopy 

Optical long-slit sp ectros copy for the UCM survey was ob- 
tained by Gallego (1995) at the 2.5-m Isaac Newton Tele- 
scope (INT) at Roque de los Muchachos Observatory, La 
Palma (Spain), and the 2.2-m and 3.5-m telescopes at Calar 
Alto (Almeria, Spain), during a total of 10 observing runs. 
Details about the instrumental setups, slit widths, spatial 
scales and dispersions achieved are given in Table 1 of 
GAL96. 

The line fluxes and equivalent widths of different 
emission lines are given in GAL96. Gas colour excesses, 
E(B — U) gas , were obtained assuming a Galactic ex- 
tinction curve and intrinsic ratios I(Ha)/I(H/3)=2.86 and 
I(H7)/I(H/3)=0.468, which are the theoretical values ex- 
pected for a low density gas w ith T C = 10 4 K in Case B re- 
combination (Osterbrock 1989). We estimate that the effect 
of differential atmospheric refraction on the Ha/H/3 ratio is 
in most cases (82% of the galaxies) below 5%. Only in four 
of the galaxies studied here the uncertainty in E(B — V) sas 
due to differential refraction is larger than 0.1 m . The obser- 



vation and analysis procedures followed by Gallego ( 1995 ) 
— slit widths, position angles, spectrum extraction — ensure 
good integrated spectroscopic information for these objects. 
Errors in the EW(Ha+[NII]) have been estimated from the 
signal-to-noise and spectral resolution data given by GAL96. 
We have assumed a lOOA interval for the continuum fit 
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Table 1. Optical-nIR colours at three disk-scale length apertures (dL=disk-scale radius, as given by Vitores et al., 1996a), and integrated 
ftf-band magnitudes. Redshifts, EW(Ha+[NII]), Lh oi E(B — V) gas and spectroscopic types have been taken from GAL96. In this table 
only the Ha luminosity data are given corrected for extinction. 

Galaxy Redshift rf L (kpc) r-J J-K K EW 1 " (A) L Ha (10 s Lq) E(B - V) gas Type 



0003+2200 


0.0245 


1.66 


1.59+0.28 


1.05+0.44 


12.83+0.35 


50+1 


0.28 


0.87 


DANS 


0013+1942 


0.0270 


2.02 


1.54+0.08 


0.95+0.10 


14.11+0.07 


142+4 


0.86 


0.28 


HUH 


0014+1748 


0.0182 


13.76 


2.07+0.11 


0.93+0.09 


11.08+0.05 


135+2 


2.95 


0.81 


SBN 


0014+1829 


0.0182 


0.88 


1.32+0.18 


0.95+0.24 


12.96+0.20 


146+2 


0.56 


1.47 


HUH 


0015+2212 


0.0199 


1.31 


1.70+0.11 


0.94+0.10 


13.21+0.07 


147+2 


1.01 


0.22 


HUH 


0017+1942 


0.0259 


2.87 


1.38+0.13 


0.84+0.11 


13.11+0.07 


181+6 


2.96 


0.36 


HUH 


0022+2049 


0.0185 


1.57 


2.15+0.12 


1.14+0.10 


11.19+0.05 


106+2 


1.83 


0.90 


HUH 


0050+2114 


0.0245 


2.39 


2.00+0.12 


0.99+0.13 


11.65+0.09 


111+1 


2.76 


0.81 


SBN 


0145+2519 


0.0409 


5.81 


1.77+0.13 


1.10+0.14 


12.11+0.10 


38+1 


2.60 


1.03 


SBN 


1255+3125 


0.0252 


1.89 


1.90+0.15 


0.88+0.23 


12.50+0.18 


74+1 


1.45 


0.41 


HUH 


1256+2823 


0.0315 


2.58 


1.56+0.11 


1.15+0.15 


12.46+0.11 


109+2 


2.82 


0.64 


SBN 


1257+2808 


0.0171 


1.11 


1.35+0.33 


1.35+0.44 


12.78+0.29 


42+1 


0.29 


1.34 


SBN 


1259+2755 


0.0240 


2.65 


1.42+0.17 


1.10+0.18 


11.89+0.13 


62+1 


1.81 


0.91 


SBN 


1259+3011 


0.0307 


2.01 


1.45+0.15 


1.28+0.19 


12.63+0.14 


34+1 


0.75 


0.68 


SBN 


1302+2853 


0.0237 


1.55 


1.73+0.15 


0.92+0.25 


12.85+0.20 


48+1 


0.43 


0.62 


DHIIH 


1304+2808 


0.0205 


2.84 


1.63+0.14 


1.26+0.19 


12.02+0.14 


33+1 


0.55 


0.11 


DANS 


1304+2818 


0.0243 


2.82 


1.43+0.11 


1.12+0.12 


12.42+0.10 


115+11 


2.31 


0.11 


SBN 


1306+2938 


0.0209 


1.64 


1.45+0.10 


1.15+0.11 


12.15+0.09 


133+4 


2.11 


0.50 


SBN 


1307+2910 


0.0187 


5.81 


1.71+0.34 


1.21+0.43 


10.37+0.28 


39+1 


2.56 


0.97 


SBN 


1308+2950 


0.0242 


7.58 


2.10+0.10 


1.17+0.14 


10.75+0.10 


59+1 


3.08 


1.38 


SBN 


1308+2958 


0.0212 


4.17 


1.71+0.06 


0.77+0.16 


12.03+0.15 


26+1 


7.03 


1.31 


SBN 


1312+2954 


0.0230 


2.97 


1.90+0.14 


1.06+0.35 


12.15+0.32 


65+3 


0.96 


1.09 


SBN 


1312+3040 


0.0210 


2.31 


1.83+0.10 


1.13+0.09 


11.69+0.07 


81+2 


1.51 


0.47 


SBN 


1428+2727 


0.0149 


1.48 


0.81+0.18 


0.84+0.22 


12.45+0.17 


218+3 


2.35 


0.15 


HUH 


1432+2645 


0.0307 


6.47 


1.71+0.10 


1.15+0.14 


11.84+0.10 


47+1 


2.09 


0.91 


SBN 


1440+2511 


0.0333 


4.40 


1.71+0.05 


1.32+0.13 


12.89+0.12 


35+1 


0.57 


1.02 


SBN 


1440+2521N 


0.0315 


2.55 


1.86+0.32 


1.31+0.43 


12.53+0.29 


104+3 


1.59 


0.77 


SBN 


1440+2521S 


0.0314 


2.06 


1.40+0.33 


1.46+0.45 


13.25+0.30 


100+5 


1.04 


0.29 


SBN 


1442+2845 


0.0110 


1.27 


1.96+0.10 


0.98+0.14 


11.68+0.10 


135+3 


0.66 


0.68 


SBN 


1443+2548 


0.0351 


3.05 


2.00+0.36 


0.61+0.44 


12.59+0.26 


76+1 


2.63 


0.73 


SBN 


1452+2754 


0.0339 


2.74 


2.40+0.37 


0.88+0.44 


12.13+0.25 


135+2 


3.09 


0.73 


SBN 


1506+1922 


0.0205 


3.00 


2.12+0.36 


1.00+0.44 


11.78+0.25 


140+6 


1.95 


0.45 


HUH 


1513+2012 


0.0369 


2.20 


1.76+0.09 


1.49+0.08 


11.87+0.07 


150+2 


6.18 


0.54 


SBN 


1557+1423 


0.0275 


1.85 


1.82+0.11 


1.01+0.09 


12.92+0.06 


54+1 


0.64 


0.37 


SBN 


1646+2725 


0.0339 


1.81 


1.72+0.23 


0.93+0.17 


15.04+0.13 


225+3 


0.50 


0.29 


DHIIH 


1647+2729 


0.0366 


3.16 


1.94+0.10 


0.99+0.10 


12.42+0.06 


59+1 


2.06 


0.89 


SBN 


1647+2950 


0.0290 


2.97 


1.88+0.34 


0.91+0.43 


11.91+0.28 


110+2 


3.76 


0.74 


SBN 


1648+2855 


0.0308 


1.95 


1.25+0.08 


1.06+0.12 


12.67+0.12 


240+15 


6.23 


0.25 


HUH 


1654+2812 


0.0348 


2.57 


1.55+0.14 


0.97+0.22 


14.93+0.18 


70+3 


0.33 


0.31 


DHIIH 


1656+2744 


0.0330 


1.04 


1.95+0.17 


1.25+0.14 


13.08+0.09 


108+1 


1.08 


0.58 


SBN 


1657+2901 


0.0317 


1.33 


1.60+0.13 


1.29+0.15 


13.37+0.12 


80+1 


0.68 


0.56 


DANS 



Note: t Equivalent width of Ha+[NII] 



range, AA cont (Gallego 1995), and a reciprocal dispersion 
of AA~3A/pixel. Thus, 



AEW = 



snrVn 



\/EW 2 + FWZI 2 



(2) 



where SNR is the signal-to-noise ratio of the continuum, N is 
the number of points used to determine the mean continuum 
flux, i.e. N=AAcont/AA~30, and FWZI is the Full Width at 
Zero Intensity. The FWZI was computed as two times the 
Full Width at Half Maximum (FWHM) of the comparison 
arc lines. Typical FWHMs are about 12.5A (for a 3" wide 
slit with the R300V+IDS+Tek3 configuration). Equivalent 
widths (EWs hereafter) of Hq+[NII] and their correspond- 
ing errors are given in Table [j]. The Ha equivalent width 
data used in this work were corrected for contamination of 
the [NII]A6548A and [NII]A6584A emission lines using the 
[NII]AA6548,6584A/Ha line ratios given by GAL96 (see also 



Gallego 1995). Finally, we consider a correction of 2 A in 
the EW(Ha) due to the Ha underlying absorption in G-K 
giants (see, e.g., Kennicutt 1983). 



4 EVOLUTIONARY SYNTHESIS MODELS 

Although the if-band luminosity can provide a very good 
estimate of the stellar mass in galaxies (Aragon-Salamanca 
et al. 1993; Chariot 1998), the contribution from red su- 
pergiants associated with recent star forming events may 
lead to the overestimate of the stellar mass when standard 
mass-luminosity relations are used. Thus, in order to es- 
timate the relative contribution of the old underlying and 
young stellar populations to the magnitudes and colours 
measured, we have developed a complete set of evolutionary 
synthesis models. These models are based on those devel- 



© 1999 RAS, MNRAS 000, 



Star formation properties of local galaxies 5 

Table 1. (cont.) 



Galaxy Redshift d L (kpc) r-J J-K K EW+ (A) L Ua (10 8 L & ) E(B — V) gas Type 



2238+2308 


0.0238 


4.44 


1.96+0.09 


0.97+0.09 


11.10+0.05 


69+1 


3.15 


1.05 


SBN 


2239+1959 


0.0242 


2.79 


1.62+0.10 


1.03+0.08 


11.52+0.04 


173+5 


5.00 


0.54 


HUH 


2250+2427 


0.0421 


6.01 


1.82+0.10 


1.20+0.08 


11.72+0.04 


175+4 


11.12 


0.71 


SBN 


2251+2352 


0.0267 


1.13 


1.53+0.11 


0.97+0.08 


13.27+0.05 


84+1 


0.97 


0.18 


DANS 


2253+2219 


0.0242 


1.45 


1.99+0.08 


1.10+0.08 


12.40+0.05 


86+1 


1.07 


0.54 


SBN 


2255+1654 


0.0388 


5.75 


2.43+0.09 


1.37+0.10 


11.56+0.05 


40+1 


1.39 


1.47 


SBN 


2255+1926 


0.0193 


2.13 


1.47+0.13 


0.96+0.11 


13.56+0.07 


37+1 


0.16 


0.37 


DHIIH 


2255+1930N 


0.0189 


2.11 


2.00+0.10 


1.07+0.09 


11.63+0.05 


97+1 


1.38 


0.70 


SBN 


2255+1930S 


0.0189 


1.12 


1.86+0.10 


0.99+0.09 


12.74+0.05 


61+1 


0.52 


0.49 


SBN 


2258+1920 


0.0220 


1.87 


2.04+0.10 


0.98+0.10 


12.35+0.06 


190+2 


1.75 


0.35 


DANS 


2300+2015 


0.0346 


2.42 


2.05+0.12 


1.10+0.10 


12.70+0.05 


159+1 


3.17 


0.33 


SBN 


2302+2053E 


0.0328 


3.09 


1.76+0.12 


1.10+0.10 


11.55+0.05 


36+1 


1.68 


1.30 


SBN 


2302+2053W 


0.0328 


2.00 


1.73+0.12 


0.97+0.11 


14.26+0.08 


260+2 


1.38 


0.46 


HUH 


2303+1856 


0.0276 


2.35 


2.10+0.13 


1.16+0.10 


11.35+0.05 


79+1 


2.41 


1.20 


SBN 


2304+1640 


0.0179 


1.00 


1.30+0.15 


0.93+0.19 


14.82+0.15 


155+2 


0.20 


0.33 


BCD 


2307+1947 


0.0271 


1.63 


2.10+0.22 


1.10+0.14 


12.41+0.08 


45+1 


0.64 


0.45 


DANS 


oqi 04-1 son 


U.UJDJ 


2 84 




i i R4-n 1 a 


±z.zuinu.uo 


UOIC-L 


1 46 


u.yu 


SBN 


2313+1841 


0.0300 


2.42 


1.99+0.12 


0.97+0.14 


13.03+0.09 


82+2 


0.72 


0.91 


SBN 


2316+2028 


0.0263 


1.42 


2.75+0.14 


1.03+0.14 


12.81+0.09 


99+4 


0.49 


0.75 


DANS 


2316+2457 


0.0277 


3.77 


2.01+0.12 


1.14+0.14 


10.41+0.08 


109+1 


10.58 


0.69 


SBN 


2316+2459 


0.0274 


4.08 


2.17+0.13 


0.94+0.14 


11.90+0.08 


72+1 


10.16 


0.57 


SBN 


2319+2234 


0.0364 


2.69 


2.46+0.12 


1.05+0.14 


12.98+0.09 


108+1 


1.91 


0.59 


SBN 


2321+2149 


0.0374 


2.74 


1.70+0.13 


0.92+0.14 


13.27+0.09 


68+1 


1.38 


0.56 


DANS 


2324+2448 


0.0123 


3.14 


2.21+0.12 


0.99+0.14 


9.60+0.08 


10+1 


1.20 


0.26 


SBN 


2327+2515N 


0.0206 


1.39 


1.20+0.12 


0.92+0.21 


13.10+0.18 


289+7 


1.83 


0.47 


HUH 


2327+2515S 


0.0206 


1.80 


1.22+0.12 


0.92+0.19 


12.82+0.15 


104+1 


1.07 


0.36 


HUH 



Note: t Equivalent width of Ha+[NII] 



oped by AH96, but use the new population synthesis models 
of Bruzual & Chariot (private communication, BC96 here- 
after), instead of the old Bruzual & Chariot (1993; BC93 
hereafter) models. From the number of ionizing photons 
supplied by the BC96 models, we have also calculated the 
contribution of the hydrogen and helium emission-lines and 
nebular continuum to the optical and nIR passbands. 

AH96 demonstrated that the properties of most of the 
UCM star-forming galaxies are better reproduced with in- 
stantaneous burst models rather that models with constant 
star formation. Therefore, we have computed the evolution 
with time of the optical-nIR colours and EW(Ha) of an 
instantaneous burst superimposed on a 15 Gyr old evolv- 
ing population. A Scalo Initial Mass Function (IMF; Scalo 
1986) with lower and upper mass cutoffs of Mi ow =0.1Mq 
and A# up =125Mq was adopted. The Cousins- R magnitudes 
given by the BC96 models have been converted into Gunn- 
r magnitudes using the relation r=_Rc+0.383— 0.083x (V — 
Rc) (Fernie 1983; Kent 1985). 

In order to match the colours predicted by the BC96 
models for a 15 Gyr old Single Stellar Population (SSP here- 
after; r — J=2.09, J — -fsf=0.85) to those measured in the 
bulges of local relaxed spiral galaxies, assuming a negligi- 
ble dust reddening (see Peletier & Balcells 1996; Fioc & 
Rocca-Volmerange 1999), we applied a small correction to 
our models 

(r - J) obs =(r - J)™ G d yr - 0.03 
(J - K)° bs =(J - if)f 5 ° G d yr + 0.06. 

In addition, the stellar mass-to-light ratio predicted by the 
model for a 15 Gyr old stellar population, 1.34 Mq/L^ q, 
was corrected to match that measured in local relaxed spiral 



galaxies, ~1Mq/L k q (see Heraudeau fc Simien 1997, and 



references therein), using M K q=3.33 ( Worthey 1994). Since 
most UCM galaxies (a bout 83 per c ent) are morphologically 
classified as Sa-Sc+ ( Vitores 1994 ), we are confident that 
the corrections applied to the models are reasonable. In any 
case, these small corrections, intended to provide a good 
agreement between the model predictions and observations 
for the underlying stellar populations of the galaxies, do not 
affect significantly any of the conclusions of this paper. 

The main parameters of our models are those inherent 
to the BC96 models (age, metallicity, IMF, etc), together 
with the strength of the current star-forming burst. The 
burst strength, b, is defined as the ratio of the mass of the 
newly formed stars to the total stellar mass of the galaxy 
(Kriiger, Fritze-v. Alvensleben & Loose 1995). We have ex- 
plored models with metallicities between 1/50 Zq and 2 Zq, 
and burst strengths in the range 1-10 -4 (in 0.04 dex steps). 
The new BC96 models (Scalo IMF) produce slightly redder 
r — J and J — K colours and fewer Lyman photons than the 
BC93 ones (Salpeter IMF), for the same burst strength and 
solar metallicity. 



5 RESULTS 

The main goal of this work is the characterization of the star 
formation activity of a representative sample of local galax- 
ies. The properties of the current star-formation events and 
the host galaxies will be studied. In particular, we are in- 
terested in linking the properties of the local star-forming 
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galaxies with those of galaxies forming stars at higher red- 
shifts. 

First, we study the completeness and representativeness 
of our sample in relation to the local star- forming galaxy 
population (see section 5.1). In section 5.2 we analyze the 
measured magnitudes and colours of the galaxies. Then, and 
in order to obtain the burst strengths, burst ages, stellar 
masses and specific star formation rates (SFR per unit mass; 
Guzman et al. 1997; Lowenthal et al. 1997), we co mpa re our 
data with evolutionary synthesis models (section 5.3). The 
comparison between data and models, and the determina- 
tion of the best-fitting set of parameters, are not straight- 
forward tasks, and some details of the method we have ap- 
plied are described in the Appendix. Finally, we discuss the 
derived burst strengths, ages, metallicities, galaxy stellar 
mas ses and s tar form ation rates derived for our sample (sec- 



tion 5.4, 



5.5 and 



5.1 Sample completeness 

Our nIR sample will suffer, first, from the intrinsic selection 
effects of the objective-prism+photographic plate technique 
used in the UCM survey. Those were di scusse d in detail in 
Zamorano et al. (1994, 1996), Vitores (|l994|) and Gallego 



(1995). Briefly, the observational procedure employed limits 



the UCM sample to local galaxies with redshift lower than 
0.045±0.005 and Ha equivalent width larger than 20A. The 
Gunn-r limiting magnitude is 16. 5 m with a brigth end cut- 
off, due to saturation of the photographic plates, placed at 
~14.2 m . However, additional selection effects may be present 
in our work due to the limited size of our nIR sample (~35 
per cent of the UCM whole sample), thus it is necessary to 
ensure that the properties of this subsample are representa- 
tive of those of the complete UCM survey. 

In Figure [l]a we compare the Gunn-r apparent magni- 
tude histogram of the whole UCM sample with that of the 
galaxies observed in the nIR. Although the apparent magni- 
tude distributions match reasonably well, the objects in the 
nIR subsample tend to be marginally brighter than those in 
the UCM complete sample. The median m r for the UCM 
complete sample is 15.5 m , while that of the nIR sample is 
15. 2 m (see Vitores et al. 1996b). This may imply a small 
deficiency of low- luminosity and/or higher redshift galax- 
ies. However, that does not seem to be a strong effect (see 
figures ^b and d) . 

From Figures |l]b, ^c and|l|e it is clear that the nIR sam- 
ple represents abo ut 35 per cent of the UCM complete sam- 
ple flGallego 1995| ) in every redshift, E(B-V) and EW(Hq) 
bin. A Kolmogorov-Smirnov test indicates that both samples 
show similar distributions in z, E(B~V) and EW(Hq), with 
probabilities of 45, 93 and 87 per cent, respectively. For the 
r-band absolute magnitudes and Ha luminosities the prob- 
abilities are 25 and 47 per cent, respectively (see Figures |l|d 
and |l|f). Thus the galaxies in the nIR sample seem to be a 
fair subsample of the UCM complete sample in their global 
properties. 

The only small difference arises when comparing in de- 
tail the spectroscopic type distributions of the star-forming 
galaxies in the nIR and UCM complete samples. There is a 
small deficiency of Hll-hke galaxies (see ahead) in the nIR 
subsample relative to the whole UCM sample. About 30 per 
cent of the UCM whole sample are Hll-like galaxies, whereas 
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Figure 1. From top to bottom, distributions of the observed 
Gunn-r magnitudes, redshifts, gas colour excesses, absolute 
Gunn-r magnitudes, EW(Ha), and Ha luminosities. The open 
histograms correspond to the complete UCM sample, and the 
grey-filled areas correspond to the nIR sample. 
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only 19 per cent are present in the nIR subsample. Conse- 
quences of such a limitation will be taken into account in 
further discussions. 



5.2 Aperture magnitudes and colours 

5.2.1 Mean colours 

Global colours (obtained inside three disk-scale lengths), 
together with the integrated Jf-band magnitudes (not cor- 
rected for extinction) and the E(B — V) gas colour excesses, 
are given m Table m. In GAL96, the galaxies in the UCM 
sample were classified in different morphological and spec- 
troscopic classes (listed in table ^). We will briefly describe 
them here (see GAL96 for details): 

SBN — Starburst Nuclei — Originally defined by Balzano 
(1983), they show high extinction values, with very low 



Table 2. Mean colours, Ha equivalent widths (A), and corre- 
sponding standard deviations of the mean for the UCM nIR sam- 
ple. 



[Nil] /Ha ratios and faint [OIII]A5007 emission. Their Ha 
luminosities are always higher than 10 8 Lq . 
DANS — Dwarf Amorphous Nuclear Starbu rst — Intro- 
duced by Salzer, MacAlpine & Boroson ( 1989 ), they show 
very similar spectroscopic properties to SBN objects, but 
with Ha luminosities lower than 5xl0 7 Lq. 
HUH — HII Hotspot — The HII Hotspot class shows (see 
GAL96) similar Ha luminosities to those measured in SBN 
galaxies but with large [OIII] A5007/H/3 ratios, that is, higher 
ionization. 

DHIIH —Dwarf HII Hotspot— This is an HUH subclass 
with identical spectroscopic properties but Ha luminosities 
lower than 5x10 7 Lq. 

BCD — Blue Compact Dwarf — Finally, the lowest lumi- 
nosity and highest ionization objects have been classified 
as Blue Compact Dwarf galaxies, showing in all cases Ha 
luminosities lower than 5x10 7 Lq. They also show large 
[OIII]A5007/H/3 and Ha/[NII]A6584 line ratios and intense 
[OII]A3727 emission. 

In our analysis, we separate the galaxies in two main 
categories: starburst disk-like (SB hereafter) and Hll-like 
galaxies (see Guzman et al. 1997; Gallego 1998). The SB-/ifce 
class includes SBN and DANS spectroscopic types, whereas 
the Yill-Uke includes HUH, DHIIH and BCD type galaxies. 

In order to determine representative mean optical-nIR 
colours for each galaxy group, we have assumed Gaussian 
probability distributions for the r — J and J — K colours 
and EW(Ha) with the centres and widths (a) given in Ta- 
ble We have weighted the data for each galaxy with their 
corresponding errors when determining the mean values. 

The Hll-like objects seem to be on average 0.2 m bluer 
in r — J and 0.1 m in J — K than the SB galaxies (see Ta- 
ble ^). Sinc e the mean colour excess of the SB population 
(E(B- V)=0.7 m ) is 0.2 m higher than that of the HII-Mie 
galaxies, these colour differences are even more significant 
when data not corrected for extinction are used: the dif- 
ferences for the un-corrected colours are 0.35 m in r — J and 
0.15 m in J-K. K-S tests performed on the SB-hke and HII- 
like objects indicate that both subsamples arise from inde- 
pendent distributions with a probability of 99.7 and 99.9 per 
cent respectively for the un-corrected r—J and J—K colours. 

Finally, whereas more than 60 per cent of the Hll-like 
objects show not corrected for extinction equivalent widths 
of Ha higher than 120 A, only 3 per cent of the SB galaxies 
do. The relatively low EW(Ha) detection limit estimated 





n r — J a J — K a 


EW(Hq) a 


Not corrected for extinction 


Total 


67 1.79 0.05 1.06 0.04 


80 7 


SB-like 


49 1.88 0.05 1.10 0.05 


60 5 


(SBN + DANS) 






mi-like (HiiH+ 18 1.54 0.09 0.96 0.30 


133 16 


DHIIH + BCD) 






Corrected for extinction 


Total 


66 1.26 0.05 0.87 0.06 


168 10 


SB-like 


49 1.28 0.06 0.89 0.09 


150 10 


Hll-like 


17 1.21 0.07 0.81 0.35 


220 30 



for the UCM survey (~20A; Gallego 1995) ensures that the 
difference in EW(Ha) between SB and HII-/ifce galaxies is 
not due to selection effects. A K-S test gives a probability of 
99.9 per cent that these samples have independent EW(Ha) 
distributions. These differences in colours and Ha equiva- 
lent widths are probably related to differences in their evolu- 
tionary properties (typical starburst age, starburst strength, 
specific star form ation rate, etc) between both galaxy types 
(see section 



5.4 



and section 5.6) 



5.2.2 Colour-colour and colour-EW(Ha) diagrams 

In Figure |^ we show colour-colour {r — J vs. J — K) and 
colour- EW(Ha) plots for the nIR sample. The offset between 
the position of the star forming galaxies [filled circles) in the 
r — J-J — K plane (Fig ure ^a) and the bulges an d disks of 
relaxed nearby spirals jPeletier fc Balcells 199fj) indicates 
the existence of ongoing star formation. Error bars in Fig- 
ure |^a represent ±lcr errors. In Figures ^a and ^3 we plot 
solar metallicity models with burst strengths, 10~ 3 , 10 -2 , 
10 _1 , and 1. In the case of Figures ||c and ^d, models with 
10 _1 burst strength and different metallicities are displayed 
(cf. section fl ). 

Figures |a and |2|c show that changes in the optical- 
nlR colours due to changes in burst strength and age are 
much more significant than those produced by changes in 
metallicity. This fact is also observed in the colour-EW(Ha) 
diagrams (Figures ^b and ^|i) , especially in the case of sub- 
solar metallicity models. It is thus clear that it is in prin- 
ciple possible to infer burst strengths and ages from these 
diagrams, but the determination of metallicities would be 
very uncertain. 



5.3 Determination of the physical properties of 
the galaxies 

For each individual galaxy we have information on its 
r—J and J — K colours and Ha equivalent width. Thus, 
each galaxy has a point associated in the r — J, J — K, 
2.5xlogEW(Ha) three-dimensional space. However, due to 
the calibration and photometric errors, the uncertainty in 
these measurements transfor ms th ese points into probability 
distributions. As in section 5.2.1, we will assume Gaussian 
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Figure 2. Colour-colour and colour-EW(Ha) diagrams. In panels a and b optical-nIR colours are plotted for the 67 galaxies of the sample. 
Solar metallicity models have been plotted using progresively thicker lines for higher burst strenght models (b=10 -3 ,10 — 2 ,10 -1 ,1). Data 
for nearby relaxed spiral galaxies have been taken from Peletier & Balcells (1996); both bulge (asterisks) and disk (triangles) colours are 
shown. Panels c, d show the optical-nIR, colours of the SB-like (dots) and the HH-tefce (stars) objects. Values not corrected for extinction 
are also shown (small circles). In panels c &: d, models with 10 -1 burst strength and different metallicities from 1/50 Zq to 2 Zq have 
been drawn (1/50 Zq, thin solid line; 1/5 Zq, dashed line; 2/5 Zq, dashed-dotted line; Zq, thick solid line; 2Zq, dotted line), drier error 
bars are also shown. 



probability distributions for the r — J, J — K colours and 
2.5xlogEW(Ha) with the centres and widths (a) given in 
Table hi. The evolutionary synthesis models that we will as- 
sociate with each galaxy probability distribution, will also 
follow different tracks in this three-dimensional space. 



The three-dimensional probability distributions (r — 
J, J — iT,2.5xlog EW(Hq)) have been reproduced using a 
Monte Carlo simulation method. A total of 10 3 data points 
were generated in order to reproduce this distribution for 
each galaxy. No significant differences were observed using 
a larger number (e.g., 10 4 ) of input particles. We estimated 
the model that better reproduces the colours and EW(Hq) 
for each of the 10 3 test particles applying a maximum like- 
lihood method. The maximum likelihood estimator used, €., 
includes two colour terms and an EW(Ha) term. Thus, 



c{t,b,z) = Yl 



2tyAC, 



■ exp 



(c n C n )* 
2AC n 2 



(3) 



where Ci , C2 and C3 are the r — J and J — K colours and 
2.5xlogEW(Ha) and Ad, AC 2 and AC 3 are their corre- 
sponding errors. The c„ coefficients are the r — J, J — K 
and 2.5xlogEW(Ha) values predicted by a given mo del. A 
similar estimator was employed by A braham et al. (199!: ) 
for a sample of intermediate-,? HDF (Williams et al. 1996 ) 
galaxies. 

Finally, we obtained the age, t, burst strength, b, and 
metallicity, Z, of the model that maximizes C for each test 
particle of the (r — J, J — i(',2.5xlogEW(HQ:)) probability 
distribution. Therefore, this procedure effectively provides 
the (t,b,Z) probability distribution for each input galaxy. 

The resulting (t,b,Z) probability distributions are in 
many cases multi-peaked. Instead of analyzing these proba- 
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bility distributions as a whole, we have studied the clustering 
pattern present in the (t,b,Z) solution space. We have used 
a single linkage hierarchical clustering method (see Murtagh 
& Heck 1987; see also Appendix) , which allows to isolate dif- 
ferent solutions in the (t,b,Z) space. We have recovered the 
three most representative solution clusters for each galaxy. 
In Table ^ we show the mean properties of those solution 
clusters with probability higher than 20 per cent. This prob- 
ability is computed as the number of test particles in a given 
cluster over the total number of test particles (10 3 ). The er- 
rors shown in Table ^ correspond to the standard deviation 
of the data for each solution cluster. In those cases where 
all the solutions within a cluster yield the same age, burst 
strength, metallicity or mass, no errors were given. 

The subsequent statistical analysis of each of the (t,b,Z) 
clusters indicates that significant correlations between t, b 
and Z are present. We have performed a principal compo- 
nent analysis (PCA hereafter; see Morrison 1976; see also 
Appendix) of the individual clusters given in Table |j. The 
orientation of the first PCA component and the contribution 
of this component to the total variance within the solution 
cluster are also given. 

After applying this procedure to the observed sample, 
only three galaxies (UCM1440+2521S, UCM1506+1922 and 
UCM1513+2012; see Figure |a and |b) show £ max <10.0. 
Note that a value £ max = 10.0 corresponds to a model where 
the differences between the observed data and the model 
predictions equal the measurement errors, assuming mean 
errors of 0.12 m in (r — J) and (J — K) and 10 per cent in 
EW(Ha) . This indicates that the range of model predictions 
covers reasonably well the observed properties of the galax- 
ies. In Figure and ^b we plot the differences between the 
colours and EW(Ha) measured and those predicted by the 
best-fit model. These differences have been calculated for 
the central values of the r — J, J — K and 2.5xlogEW(Ha) 
probability distributions. 

Figure ^ shows that, in some cases, the models pre- 
dict bluer J — K colours than those measured. At first 
sight, these discrepancies could be explained if we were un- 
derestimating the extinction correction factors applied to 
these objects. However, since the extinction correction af- 
fects r — J more than J — K, applying a higher extinc- 
tion correction would destroy the good agreement between 
the observed and model r — J colours. We have quantified 
the effect of a change in the correction for extinction as- 
sumed for our galaxies by comparing data de-reddened us- 
ing E(B — V)continuum=-E(-B — V) gas with our models. Fig- 
ure ^ shows that data corrected using the relation given 
by Calzetti et al. (1996) fit the models better than data 
corrected assuming E(B — l / ) C ontinuum=i5(S — V) gas . In ad- 
dition, none of the three galaxies given above show higher 
£max values after using E(B - V) C ontinuum=-B(B - ^) gas . 
Thus, we are confident that the extinction correction applied 
provides a reasonable fit to the models, and the discrepan- 
cies in J — K between the data and the models are probably 
due to inherent uncertainties in the modelling of the nIR 
continuum by the Bruzual and Chariot code. 

The cluster analysis performed indicates that the clus- 
tering in the solution space is basically produced by the 
discretization in metallicity of the models. Fortunately, in 
many of the objects (~30 per cent; see Table ^), only one 
(t,b,Z) solution cluster is able to reproduce the observables. 
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Figure 3. a) & b) Differences between the r — J , J — K and 
EW(Hq) measured values and the best-fitting model values, us- 
ing E(B-V) continuum =0.UxE(B-V) sas . The input r-J, J-K 
and EW(Ha) data correspond to the central values of the respec- 
tive probability distributions. Those galaxies with £ max <10 are 
shown as grey dots, c) Differences in the r — J and J — K colours 
assuming E{B - V) CO ntinuum=-E ; (B - V r ) gaa . 



About 33 per cent need two solutions and three solutions 
are needed for the remaining 37 per cent. The goodness 
of this comparison method, given by the number and size 
of statistically significant (t,b,Z) solution clusters, basically 
depends on the particular position of the object in the 
(r — J, J — if,2.5xlogEW(Ha)) space and on its measure- 
ment errors. Fortunately, in those cases where the (t,b,Z) 
probability distribution is multi-valuated, the different so- 
lution clusters give similar burst strengths and total stellar 
masses. This is another manifestation of the fact that, as 
our data is not very sensitive to 



we saw in section 



5.2.2 



metallicity, and we will not attempt to derive it. Nonethe- 
less, the burst ages are affected somewhat by small changes 
in metallicity, and frequently show wider distributions than 
the burst strenghts. 
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Table 4. Mean properties and standard deviations for the nIR 
sample. Dwarfs includes the DHIIH and BCD spectroscopic type 
galaxies. 



t a 



logb 



logZ cr logM 
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with E(B - V) 



continuum ■ 



=0A4xE(B - V) g 



Total 


67 


5.5 


0.4 


-1.72 


0.07 


-0.5 


0.1 


10.69 


0.06 


SBN 


11 


5.8 


0.5 


-1.69 


0.10 


-0.5 


0.1 


10.90 


0.06 


DANS 


8 


5.7 


1.6 


-1.91 


0.23 
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0.3 


10.57 
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HUH 


13 


1.1 


0.6 


-1.72 


0.15 


-0.3 


0.2 
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0.10 


Dwarfs 


5 


6.6 


1.7 


-1.62 


0.20 


-0.7 


0.4 


9.93 


0.15 


SB 


49 


5.8 


0.5 


-1.73 


0.09 


-0.6 


0.1 


10.85 


0.05 


HII 


18 


1.7 


0.7 


-1.69 


0.10 


-0.4 


0.2 


10.29 


0.10 



Total 

SB 

HII 



with E(B - V) 



continuum 



E(B - V) 



gas 



67 
49 
18 



11.5 
12.5 
8.3 



0.6 
0.6 
1.0 



-0.77 
-0.69 
-1.01 



0.07 -1.2 
0.08 -1.4 
0.11 -0.7 



0.1 10.64 0.05 
0.1 10.77 0.05 
0.2 10.26 0.11 
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Finally, the PCA performed on each solution clus- 
ter suggests that the best-axis, given by the vector 
(wt,wiog(b):Miog(Z))=(wx,W!/,Mz) shown in Table [| is com- 
monly placed in the u x -u y (t-b) plane and obeys u x ~u y . This 
implies that age and burst strength are in many cases degen- 
erated and, therefore, the properties of an individual object 
can be reproduced both with a young, low burst strength or 
an old, high burst strength model, within the ranges given 
in Table 0. 
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5.4 Burst strengths and ages 

In Table ^| we give mean burst strengths and ages for the in- 
dividual solutions with probability higher than 20 per cent. 
Errors given are the standard deviations of the data points in 
each solution. For the stellar masses, the error related with 
the uncertainty in the A"-band absolute magnitude determi- 
nation is also given (in parenthesis) . Using these probability 
distributions we have derived the burst strength, age, mass 
and metallicity frequency histograms for the whole sample as 
well as for the SIS- like and HII- like galaxies (see Figures ^a- 
d). The number of points in the y axis of these figures cor- 
responds to the number of Monte Carlo test particles with 
a given burst strength, age, mass or metallicity within the 
accepted high-probability solutions. 

This analysis yields a typical burst strength of 
2xl0~ 2 with approximately 90 per cent of the sam- 
ple having burst strengths between 10 -3 and 10 — l . 
Only seven objects in the sample show burst strengths 
higher than 10~ with a probability larger than 50 per 
cent, UCM0003+2200, UCM0145+2519, UCM1257+2808, 
UCM1259+3011, UCM1308+2958, UCM1432+2645 and 
UCM1440+2511. 

Although the properties of the local star- forming galax- 
ies seem to be well reproduced with an episodic star forma- 
tion history (see also AH96), some of these objects may have 
evolved under more constant star formation rates (Glaze- 
brook et al. 1999; Coziol 1996). In those objects the in- 
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Figure 4. Frequency histograms of the derived physical prop- 
erties; a) Burst strength, b) age, c) X-band stellar mass, and 
d) specific SFR (SFR per unit mass). Upper panels show the 
histograms for the whole sample (open histograms) and for the 
SB-Kfce objects (light- grey filled histograms). Lower-panels are the 
frequency histograms of the HII-Kfce galaxies (dark-grey filled his- 
tograms) . 



stantaneous burst assumption could yield very high burst 
strengths. 

The burst strength histograms shown in Figure ^a give 
typically larger burst strengths for the Hll-like objects, es- 
pecially for the DHIIH and BCD type galaxies, than for the 
SB-like (see also Table ^). This segregation in burst strength 
is probably related to the difference in mean EWHJa) 
pointed out in section 5.2.1 (see Table ^). In Table ^ we 
also show the burst strengths and ages derived under the 
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Table 3. Best fit model results for the nIR sample. 
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1440+2511 


96.7 


12.60zb0.03 


-0.882+0.030 


-1.70 


10.78 


;o.o5) 


(+0 


707 


+0.707,-0.004) 


47.9 


1440+2521N 


42.8 


3.58±0.74 


-1.724+0.502 


+0.40 


10.80+0.25 


;o.i2) 


(+0 


707 


+0.707,+0.000) 


64.6 




33.2 


4.95±1.20 


-1.871+0.458 


-0.33+0.29 


10.86+0.08 


;o.i2) 


(+0 


704 


+0.664,-0.250) 


64.6 




24.0 


8.09±3.17 


-1.660+0.661 


-1.70 


10.90+0.03 


;o.i2) 


(+0 


707 


+0.707,+0.000) 


65.8 


1440+2521S 


61.6 


3.88±0.66 


-1.788+0.352 


+0.36+0.11 


10.52+0.15 


;o.i2) 


(+0 


720 


+0.671,-0.178) 


62.2 




31.8 


9.89±2.78 


-1.295+0.621 


-1.70 


10.57+0.03 


;o.i2) 


(+0 


707 


+0.707,+0.000) 


65.6 


1442+2845 


59.0 


4.55±0.63 


-1.986+0.193 


-0.41+0.31 


10.32+0.01 


;o.o4) 


(+0 


738 


+0.532,-0.415) 


60.6 




20.6 


7.47±1.42 


-1.705+0.302 


-1.70 


10.32 


;o.o4) 


(+0 


707 


+0.707,+0.000) 


66.4 




20.4 


2.70±0.42 


-2.149+0.130 


+0.40 


10.32+0.01 


;o.o4) 


(+0 


707 


+0.707,+0.000) 


63.2 
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Table 3. (cont.) Best fit model results for the nIR sample. 



Galaxy Prob. Age log b log (Z/Z & ) log (M/M & ) PCA Variance 

(%) (Myr) (b=My 0ung /M total ) K,«logb,«logz) (%) 



1443+2548 


62 


3 


5 


00+1 


47 


-1 


767+0 


572 


-0 


40+0 


26 


10 


89+0 


10 





10) 


(+0 


702 


+0.710,+0 


061) 


62 


4 




21 





4 


89+3 


56 


-2 


275+0 


706 


-1 


70 




10 


98+0 


03 





10) 


(+0 


707 


+0.707,+0 


000) 


65 


7 


1452+2754 


52 


3 


3 


13+1 


39 


-2 


398+0 


462 


-0 


57+0 


24 


11 


13+0 


04 





10) 


(+0 


673 


+0.704,+0 


226) 


65 


8 




37 


8 


2 


52+2 


14 


-2 


573+0 


367 


-1 


70 




11 


14+0 


01 





10) 


(+0 


707 


+0.707,+0 


000) 


65 


5 


1506+1922 


46 


6 


2 


81+2 


98 


-2 


690+0 


518 


-1 


70 




10 


82+0 


01 





10) 


(+0 


707 


+0.707,+0 


000) 


66 


1 




33 


1 


4 


05+1 


68 


-2 


252+0 


559 


-0 


44+0 


29 


10 


78+0 


10 





10) 


(+0 


698 


+0.708,+0 


107) 


63 


4 




20 


3 


3 


55+1 


36 


-1 


675+0 


894 


+0 


40 




10 


46+0 


56 





10) 


(+0 


707 


+0.707,+0 


000) 


61 


6 


1513+2012 


90 


3 


2 


75+0 


59 


-2 


137+0 


117 


+0 


31+0 


17 


11 


29 




'S> 


03) 


(+0 


644 


+0.508,-0 


572) 


76 


5 


1557+1423 


42 


5 


11 


71+0 


41 


-1 


450+0 


118 


-1 


70 




10 


61+0 


01 





02) 


(+0 


707 


+0.707,+0 


000) 


66 


2 




34 


6 


4 


39+0 


56 


-2 


167+0 


162 


+0 


32+0 


16 


10 


59+0 


03 





02) 


(+0 


695 


+0.499,-0 


518) 


67 


9 




22 


9 


7 


16+0 


54 


-1 


824+0 


165 


-0 


64+0 


12 


10 


57+0 


03 





02) 


(+0 


651 


+0.530,-0 


544) 


77 


8 


1646+2725 


35 


2 


2 


89+1 


43 


-2 


231+0 


329 


-0 


60+0 


14 


9 


93+0 


02 





05) 


(+0 


688 


+0.693,+0 


214) 


64 


3 




34 


2 


3 


00+2 


20 


-2 


248+0 


365 


-1 


70 




9 


93+0 


01 





05) 


(+0 


707 


+0.707,+0 


000) 


66 


3 




30 


6 


1 


78+0 


91 


-2 


179+0 


187 


+0 


28+0 


18 


9 


94+0 


01 





05) 


(+0 


720 


+0.681,-0 


133) 


59 


1 


1647+2729 


66 


4 


4 


66+0 


59 


-1 


090+0 


668 


+0 


31+0 


17 


10 


53+0 


51 





02) 


(+0 


514 


+0.764,+0 


390) 


56 


6 


1647+2950 


38 


2 


4 


87+1 


29 


-1 


769+0 


515 


-0 


57+0 


15 


11 


02+0 


08 





11) 


(+0 


635 


+0.679,+0 


369) 


68 


1 




36 


6 


3 


45+1 


05 


-1 


881+0 


385 


+0 


19+0 


20 


11 


05+0 


06 





11) 


(+0 


688 


+0.584,-0 


431) 


66 


5 




25 


2 


6 


21+3 


25 


-1 


848+0 


663 


-1 


70 




11 


08+0 


01 





11) 


(+0 


707 


+0.707,+0 


000) 


66 


2 


1648+2855 


92 


2 


2 


71+0 


47 


-1 


777+0 


085 


+0 


35+0 


14 


10 


78+0 


01 





05) 


(+0 


623 


+0.518,-0 


586) 


78 


7 


1654+2812 


(50 


9 


4 


30+0 


62 


-1 


756+0 


242 


+0 


34+0 


11 


9 


92+0 


11 





07) 


(+0 


697 


+0.528,-0 


486) 


67 


8 




20 


2 


6 


43+0 


47 


-1 


577+0 


211 


-0 


53+0 


15 


9 


91+0 


07 





07) 


(+0 


692 


+0.447,-0 


567) 


68 


1 


1656+2744 


55 


6 


4 


45+1 


13 


-2 


164+0 


295 


-0 


30+0 


26 


10 


72+0 


01 





04) 


(+0 


666 


+0.628,-0 


403) 


71 


4 




22 


2 


8 


15+1 


33 


-1 


773+0 


283 


-1 


70 




10 


72 




v° 


04) 


(+0 


707 


+0.707,+0 


000) 


66 


3 




22 


2 


2 


82+0 


94 


-2 


174+0 


389 


+0 


40 




10 


70+0 


14 





04) 


(+0 


707 


+0.707,+0 


000) 


63 


1 


1657+2901 


84 


3 


3 


94+0 


40 


-1 


760+0 


190 


+0 


38+0 


09 


10 


52+0 


07 





05) 


(+0 


703 


+0.526,-0 


479) 


67 


3 


2238+2308 


63 


8 


5 


67+0 


21 


-1 


593+0 


106 


-0 


65+0 


11 


11 


22+0 


02 





02) 


(+0 


749 


+0.169,-0 


640) 


54 


6 




2.3 


7 


3 


89+0 


73 


-1 


664+0 


142 


+0 


23+0 


20 


11 


22+0 


01 





02) 


(+0 


615 


+0.527,-0 


586) 


86 


2 


2239+1959 


68 


8 


3 


19+0 


65 


-1 


895+0 


130 


+0 


25+0 


19 


11 


11+0 


01 





02) 


(+0 


606 


+0.556,-0 


569) 


88 


3 


2250+2427 


87 


4 


2 


27+0 


66 


-2 


056+0 


131 


+0 


32+0 


16 


11 


49+0 


01 





02) 


(+0 


714 


+0.466,-0 


523) 


63 


2 


2251+2352 


51 


5 


4 


27+0 


40 


-1 


943+0 


113 


+0 


38+0 


09 


10 


39+0 


04 





02) 


(+0 


683 


+0.533,-0 


499) 


70 


9 




35 


4 


11 


40+0 


48 


-1 


304+0 


139 


-1 


70 




10 


43+0 


01 





02) 


(+0 


707 


+0.707,+0 


000) 


66 


2 


2253+2219 


48 





4 


65+0 


69 


-2 


152+0 


149 


-0 


21+0 


30 


10 


72+0 


01 





02) 


(+0 


643 


+0.536,-0 


546) 


80 


3 




31 


3 


8 


49+0 


73 


-1 


794+0 


167 


-1 


70 




10 


72 




;° 


02) 


(+0 


707 


+0.707,+0 


000) 


66 


6 




20 


7 


3 


46+0 


95 


-1 


945+0 


757 


+0 


40 




10 


50+0 


51 





02) 


(+0 


707 


+0.707,+0 


000) 


65 


6 


2255+1654 


63 





3 


79+0 


67 


-1 


963+0 


116 


+0 


21+0 


20 


11 


51+0 


01 





02) 


(+0 


649 


+0.475,-0 


594) 


77 


6 




30 


3 


9 


26+0 


22 


-1 


401+0 


062 


-1 


70 




11 


52 




v° 


02) 


(+0 


707 


+0.707,+0 


000) 


66 


6 


2255+1926 


73 


2 


13 


08+0 


34 


-1 


135+0 


133 


-1 


70 




10 


02+0 


01 





03) 


(+0 


707 


+0.707,+0 


000) 


65 


9 


2255+1930N 


65 


5 


3 


53+0 


74 


-2 


102+0 


152 


+0 


18+0 


20 


10 


86+0 


01 





02) 


(+0 


643 


+0.530,-0 


552) 


79 


6 




20 


1 


5 


02+0 


55 


-2 


014+0 


225 


-0 


54+0 


15 


10 


85+0 


01 





02) 


(+0 


679 


+0.709.+0 


192) 


65 


5 


2255+1930S 


19 


9 


4 


15+0 


39 


-1 


955+0 


149 


+0 


40 




10 


38+0 


06 





02) 


(+0 


707 


+0.707,+0 


000) 


66 







28 


7 


6 


15+0 


73 


-1 


853+0 


196 


-0 


42+0 


30 


10 


39+0 


03 





02) 


(+0 


661 


+0.500,-0 


559) 


75 


8 




21 


4 


10 


97+0 


57 


-1 


428+0 


154 


-1 


70 




10 


42+0 


01 





02) 


(+0 


707 


+0.707,+0 


000) 


66 


4 


2258+1920 


67 


8 


1 


32+0 


92 


-2 


692+0 


102 


-1 


70 




10 


65 




;° 


02) 


(+0 


707 


+0.707,+0 


000) 


61 


1 




24 





1 


52+1 


12 


-2 


576+0 


115 


-0 


36+0 


25 


10 


65 




v° 


02) 


(+0 


643 


+0.708.+0 


293) 


60 


6 


2300+2015 


77 


3 


1 


27+0 


58 


-2 


713+0 


062 


-1 


70 




10 


90 




[0 


02) 


(+0 


707 


+0.707,+0 


000) 


56 


5 


2302+2053E 


79 


5 


4 


34+0 


47 


-1 


349+0 


123 


+0 


35+0 


14 


11 


23+0 


06 





02) 


(+0 


613 


+0.532,-0 


584) 


86 


7 


2302+2053W 


19 


4 


2 


38+0 


95 


-2 


239+0 


142 


-0 


38+0 


27 


10 


23+0 


01 





03) 


(+0 


734 


+0.605,-0 


308) 


57 


8 




26 


3 


1 


26+0 


63 


-2 


117+0 


096 


+0 


40 




10 


23 




;o 


03) 


(+0 


707 


+0.707,+0 


000) 


58 


3 




24 


3 


3 


65+1 


20 


-2 


077+0 


253 


-1 


70 




10 


21+0 


01 





03) 


(+0 


707 


+0.707,+0 


000) 


65 


2 


2303+1856 


55 


5 


3 


47+0 


72 


-1 


856+0 


164 


+0 


19+0 


20 


11 


27+0 


01 





02) 


(+0 


608 


+0.554,-0 


569) 


88 


7 




23 


1 


8 


40+0 


55 


-1 


334+0 


155 


-1 


70 




11 


28 




r S> 


02) 


(+0 


707 


+0.707,+0 


000) 


66 


1 




21 


4 


5 


09+0 


10 


-1 


721+0 


170 


-0 


55+0 


15 


11 


26+0 


01 





02) 


(+0 


710 


+0.691,-0 


136) 


61 


9 




1 7 


a 
o 




Q7-t-0 


uu 


_ I 




-LOU 


n 
u 


ooxu 


ox 


Q 


43-1-0 


ft'"-! 



u 


UU J 


( -\-C\ 








7R 


r. 
u 




31 


7 


8 


90+0 


51 


-1 


198+0 


175 


-1 


70 




9 


45+0 


01 





06) 


(+0 


707 


+0.707,+0 


000) 


66 


1 




20 


7 


3 


17+0 


17 


-1 


710+0 


153 


+0 


40 




9 


44+0 


02 





06) 


(+0 


707 


+0.707,+0 


000) 


62 


3 


2307+1947 


50 





4 


21+1 


50 


-2 


355+0 


516 


+0 


28+0 


18 


10 


73+0 


18 





03) 


(+0 


680 


+0.704,+0 


204) 


65 


7 




33 


2 


8 


59+4 


48 


-2 


173+0 


802 


-1 


70 




10 


81+0 


01 





03) 


(+0 


707 


+0.707,+0 


000) 


66 


5 


2310+1800 


45 


2 


4 


46+1 


14 


-2 


230+0 


271 


-0 


34+0 


29 


11 


15+0 


01 





03) 


(+0 


714 


+0.673,-0 


193) 


63 


2 




28 


2 


3 


10+1 


22 


-2 


072+0 


678 


+0 


40 




10 


99+0 


36 





03) 


(+0 


707 


+0.707,+0 


000) 


64 


4 




26 


6 


6 


98+2 


85 


-2 


032+0 


520 


-1 


70 




11 


15 




v° 


03) 


(+0 


707 


+0.707,+0 


000) 


66 


1 


2313+1841 


70 


2 


5 


03+0 


65 


-1 


777+0 


184 


-0 


43+0 


29 


10 


65+0 


02 





04) 


(+0 


675 


+0.510,-0 


534) 


72 


5 


2316+2028 


99 


8 





92+0 


36 


-2 


731+0 


026 


-1 


70 




10 


65 




r s> 


04) 


(+0 


707 


+0.707,+0 


002) 


39 





2316+2457 


43 





4 


69+0 


65 


-2 


069+0 


233 


-0 


62+0 


13 


11 


63+0 


01 





04) 


(+0 


703 


+0.709,+0 


059) 


65 


1 
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Table 3. (cont.) Best fit model results for the nIR sample. 



Galaxy 


Prob. 


Age 


log b 


log (Z/Z & ) 


log (M/M Q ) 


PCA 


Variance 




(%) 


(Myr) 


(6=M youn g / M total ) 






(«t,Mlogb,«lo g z) 


(96) 




35.5 


o ncxn on 

2.yb±U.oy 


O 1 C7 + n OQO 


i n o a _Ln on 
-hU.z4±U.zU 


1 1 cQj.n m /'n 
11.0o±U.Ul(U.U4) 


(+u.b /y,+u.oyy,— U.424J 


68.9 




Zl.D 


fi f;n_!_ 1 qi 
O.OUztl.ol 


— l.yU f ±U.o4o 


1 7H 


1 1 cka ( n n/i ~\ 
11.04 (U.U4J 


/ i n 7H7 i n '7n r 7 i n nnn\ 
^-hU. ( U ( ,+U. / U ( ,H~U.UUU ) 


00.4 




71.6 


a nn j_i oo 
4.yU±l.zo 


— l.zob±U. /y4 


i n Qi;j.n i o 

+U.oo±U.lz 


in Q7+.n t.c\( c\ n/i*\ 
lU.o / ±U.oU(U.U4) 


(+U.olO,+U.OOZ,+U.44o ) 


73.4 




si 1 


i Q74-0 ^7 


— i.OOOiU.UOU 


1 70 




(-1-0 70,7 4-0 707 4-0 OOO 1 ! 


^9 


2321+2149 


61.0 


4.51+0.57 


-1.617+0.219 


+0.33+0.15 


10.62+0.15(0.04) 


(+0.715, +0.488,-0.500) 


64.7 




30.6 


6.40+0.36 


-1.553+0.226 


-0.56+0.15 


10.65+0.09(0.04) 


(+0.733, +0.152,-0.663) 


57.7 


2324+2448 


57.3 


6.96+2.30 


-1.717+1.105 


+0.33+0.16 


10.77+0.43(0.04) 


(+0.648, +0.692,+0.317) 


67.3 




25.0 


13.51+2.61 


-2.229+0.446 


-1.70 


11.24 (0.04) 


(+0.707,+0.707,+0.000) 


66.3 


2327+2515N 


35.7 


5.89+0.99 


-1.391+0.209 


-1.70 


10.25+0.01(0.07) 


(+0.707,+0.707,+0.000) 


66.1 




33.1 


2.28+0.30 


-1.602+0.125 


+0.40 


10.24+0.01(0.07) 


(+0.707,+0.707,+0.000) 


60.6 




31.2 


3.94+0.43 


-1.548+0.179 


-0.42+0.29 


10.22+0.02(0.07) 


(+0.741, +0.473,-0.476) 


58.6 


2327+2515S 


81.9 


4.25+0.46 


-1.253+0.370 


0.37+0.11 


10.15+0.29(0.06) 


(+0.711, +0.494,-0.501) 


65.5 



unrealistic assumption that the continuum extinction is a s 
high as that measured for the ionized gas (see section |5,3| ) . 

The distribution of the burst ages is shown in Figure Mb. 
Since the probability of detection increases with EW(Ha) 
in objective-prism surveys (see Garcfa-Dabo et al. 1999), 
and the EW(Ha) continuously decreases with the the burst 
age, the number of objects detected with old burst ages is 
expected to be lower than with young ages, as observed. 
This behaviour is observed at ages older than 4 Myr, both 
for the SB and Hll-like galaxies. 

However, one would expect a reasonably fiat distribu- 
tion in the number of objects with young ages if the sample 
selection depended only on the Ha equivalent width. But 
other factors such as the Ha flux and continuum luminos- 
ity play an important role (see Garcfa-Dabo et al. 1999). 
Moreover, in our models we have estimated the Ha lumi- 
nosity (I/Ha) from the number of Lyman continuum photons 
(Brocklehurst 1971) assuming that no ionizing photons es- 
cape from the galaxies. If some Lyman photons escape, the 
predicted Ha luminosity would be lower and the derived ages 
could be significantly younger. Recent studies estimate the 
fraction of Lyman photons escaping from starbur st galax- 
ies to be about 3 per cent (Leitherer et al. 1995). Bland- 
Hawthorn & Maloney (1997) estimated this quantity to be 
about 5 per cent for the Milky Way. Another feasible expla- 
nation could be that a significant fraction of these Lyman 
photons is absorbed by dust within the ionized gas (see, 
e.g., Armand et al. 1996). Both mechanisms would produce 
lower Ha equivalent widths than those predicted by the stan- 
dard super-ionizing models, and could explain the paucity of 
young star- forming bursts in Figure Qd. In this figure (dotted 
line) we also show the age distribution obtained assuming 
that 25 per cent of the Lyman continuum photons are miss- 
ing. This distribution yields a larger number of objects at 
ages younger than 3 Myr, and a very steep decay at ages 
older than 4-5 Myr. 



Finally, Bernasconi & Maeder (1996) have argued that, 



during the first 2-3 Myr in the main-sequence, stars more 
massive than 40 Mq , are still accreting mass embedded in 
the molecular cloud, and do not contribute to the ionizing 
radiation. Therefore, due to this reduction in the number of 
Lyman photons, the predicted Ha equivalent widths below 
2-3 Myr will be significantly lower and the ages deduced for 
the bursts should be younger. 



5.5 Total stellar masses 



In order to determine the total galaxy stellar mass we have 
assumed that the burst strengths and mass-to-light ratios 
derived from our models at three disk scale-length apertures 
are representative of the galaxy global properties. Thus, us- 
ing these if-band mass-to-light ratios and the total if-band 
absolute magnitudes we have obtained stellar masses for the 
whole sample. 

The inferred galaxy stellar masses derived depend, in 
principle, on four quantities: the galaxy i\-band absolute 
magnitude, burst strength, and the mass-to- light ratios of 
the burst and the old underlying population. Since the de- 
rived burst strengths are very low (~10 -2 ), the total mass- 
to-light ratios are dominated by the old stellar component. 
In fact, the ratio of the if -band luminosity of the old and 
young stellar populations is ~20 for i=4Myr, 4 for i=8Myr 
and 7 for t=15Myr (for Z=Zq). Moreover, the absolute age 
of the old stellar component has a very small effect on the 
if-band mass-to-light ratio: there is only a 0.1 dex differ- 
ence between 10 Gyr and 15 Gyr for solar metallicity, and 
the difference is even lower in the case of sub-solar metal- 
licity models. In Figure ^ we show that the derived JT-band 
stellar mass-to- light ratios span a very narrow range. Al- 
though statistically the SB- and HII-/ifce mass-to-light ratio 
distributions are different with a probability of 95.3 per cent 
(from a K-S test), the difference in absolute value is only mi- 
nor: the median mass-to-light ratios are 0.93Mq/L k q for 
the whole sample, and 0.93Mq/L x q and 0.91Mq/L x q, 
for the SB and HB-Zifce objects respectively. Consequently, 
the derived mass values mainly depend on the if-band abso- 
lute magnitude. In the upper-panel of Figure [| we also show 
t he ra nge in If -band mass-to-light ratios given by Worthey 
( 1994 ) for 12 Gyr old modeled ellipticals. Thus, we can con- 
clude that the if -band luminosity is a very good tracer of 
the stellar mass for both old stellar populations and local 
star-forming objects. 

The distribution of stellar masses is shown in Figure ^c. 
This frequency histogram indicates that a typical star form- 
ing galaxy in our Local Universe has a stellar mass of about 
5x10 10 Mq. This value is somewhat lower than the stellar 
mass expected for a local L* galaxy. Assuming M* K =— 25.1 
(for Ho=50kms _1 Mpc -1 ; Mobasher, Sharpies and Ellis, 
1993) and a if -band mass-to-light ratio of 1 Mq /L k q 
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Figure 5. _R"-band mass-to-light ratio frequency histogram. The 
upper panel shows the distribution for the whole sample (open 
histogram) and SB-like objects (light-grey filled histogram). The 
lower-panel shows the frequency histogram of the HII-Kfce galaxies 
(dark-grey filled histogram). Mass-to-light ratios for model ellip- 
ticals with different metallicities are taken from worthey (1994, 
w94). 



Table 5. Mean and median log(Ljj Q /SFR) ratios obtained for 
different metallicities. 



Metallicity 



(L Hq /SFR) 



z 


Mean 


Median 


Std.dev. 


1/50 Z Q 


40.19 


40.09 


0.58 


1/5 Z 


40.38 


40.35 


0.36 


2/5 Z Q 


40.31 


40.27 


0.29 


z 


40.24 


40.29 


0.31 


2Z 


40.19 


40.21 


0.34 


All Z 


40.23 


40.23 


0.44 



(Heraudeau & Simien 1997), the stellar mass inferred for an 
L* galaxy is about 2x10 Mq. Thus, star-forming galax- 
ies in the local universe are typically a factor 4 less massive 
than L* galaxies. 

In addition, a clear offset between the stellar mass 
histograms of the SB and HII-Hke objects is seen in Fig- 
ure [4c. The distributions of their stellar masses are centred 
at 7x 10 10 and 2x 10 10 M respectively (Figure ||c) . This dif- 
ference is even more significant, about ldex, when DHIIH 
and BCD spectroscopic types (Dwarfs in Table ^) and SBN 
galaxies are compared. A K-S test analysis of the SB-like and 
Hll-hke objects indicates that these two samples come from 
different age, burst strength and stellar mass distributions 
with probabilities 98.8, 77.1 and 99.9 per cent respectively. 

In Table ^ we also present the mean properties that 
would be obtained using E(B — V) C ontinuum = E(B — V) gss . 
In this case, although we obtain important differences in age 
and burst strength, very similar stellar masses are derived 
since the stellar mass depends mainly on the Tf-band mag- 
nitude, only weakly affected by extinction. 



5.6 Star formation rates 

Since the star formation activity in the UCM galaxies is 
better described as a succession of episodic star forma- 



Table 6. Mean SFR and specific SFR for the sample, together 
with the corresponding standard deviations of the mean. Dwarfs 
includes DHIIH and BCD spectroscopic type galaxies. 





n 


1^,,-/QFD \ 

iOg(br ix) 


o~ 


log(br H/MJ <7 






(M Q yr~ 




(lO-Hyr- 1 ) 


Total 


67 


1.52 


0.05 


1.81 0.04 


SBN 


41 


1.64 


0.06 


1.73 0.05 


DANS 


8 


1.21 


0.09 


1.66 0.09 


HUH 


13 


1.60 


0.08 


2.16 0.07 


Dwarfs 


5 


0.82 


0.09 


1.88 0.14 


SB 


49 


1.57 


0.05 


1.72 0.04 


HII 


18 


1.38 


0.10 


2.08 0.07 



tion events rather than continuous star formation (see also 
AH96), the current star formation rate (SFR) is not a mean- 
ingful quantity: the latest star formation event might have 
finished in many of the galaxies, and their current SFR 
would be zero. However, these galaxies have substantial Ha 
luminosities, and it is accepted that the Ha luminosity is a 
good measurement of the current SFR. In AH96 we showed 
that this is true, in a statistical sense, for a population of 
galaxies undergoing a series of star formation events, and we 
defined an 'effective' present-day SFR which coincides with 
the SFR we would derive if the galaxies were forming stars 
at a constant rate, producing the same mass in new stars as 
the ensemble of all the star-formation episodes (see AH96 
for details). Here we will follow the same approach. 

When estimating star formation rates (SFRs) in AH96 
we used the BC93 models and a Salpeter IMF. In the present 
work, we have used the updated BC96 models with a Scalo 
IMF. Since the number of Lyman photons (NLy) predicted 
by the old models is about 0.94 dex higher than that pre- 
dicted by the new ones (for solar metallicity and ages lower 
than 16Myr), we need to re-compute the relation between 
the Ha luminosity, I/Ha, and star formation rate. In addi- 
tion, we will investigate the change produced in this relation 
using different metallicity models. 

In order to compute the Lhq/SFR ratio, we have used a 
very similar procedure to that employed by AH96: we sim- 
ulated a population of galaxies undergoing random bursts 
of star-formation and computed their total Ha luminosi- 
ties and the mass in newly-formed stars. However, instead 
of considering a uniform age and burst strength probabil- 
ity distribution, we have considered the burst strength, age 
and metallicity distributions for our galaxy sample. We used 
67x1000 points in order to reproduce this distribution in 
our Monte Carlo simulations. The SFR was computed as 
the ratio between the stellar mass produced in the burst, 
i.e. 6xM, and the maximum age for which we could have 
detected the galaxy i n the UCM sa mple, that is, the time 
while EW(Ha)>20A jGallego 1995| ). The L Ha /SFR ratios 
obtained are shown in Figure |^ for different metallicities 
and for the total (t,b,Z) distribution. The mean, median, 
and standard deviation values are given in Table U. 

Since the changes in this ratio for different metallicity 
models are quite small, we have adopted the median value 
of the whole distribution in order to determine the 'effec- 
tive' SFR of the galaxies from our sample. The difference 
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Figure 6. log (Lh q /SFR) frequency histograms for different 
metallicities and for the whole distribution. Thick marks give the 
position of the distribution median values. 
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Figure 7. Stellar mass (in Mq) vs. specific SFR (in 10" 



a) Dynamical masses from Telles (1995), converted to stellar 
masses using a correction factor of 0.6 dex. Starburst galaxies 
from Calzetti (1997b) are also plotted. Data for spiral galaxies 
have been taken from Kennicutt (1983). b) SBN, DANS, HII, 
DHIIH and BCD UCM galaxies as classified by Gallego et al. 
(1996) are plotted. Straight lines define the limits in the Ha lu- 
minosity function of the UCM survey. 



between the value adopted here and that of AH96 is about 
1 dex, which is very close to the difference in the number of 
Lyman photons predicted by the BC93 and BC96 models, 
as expected. 

Therefore, we have evaluated the current 'effective' SFR 
using the expression 



SFR: 



M yr" 



(4) 



1.7 x 10 40 ergs" 
This expression assumes that every Lyman photon emitted 



effectively ionizes the surrounding gas. If, however, as is sug- 
gested in Section 5.4, we consider a fraction of non-ionizing 



Lyman photons of 25 per cent, the star formation rates com- 
puted should be 0.1 dex higher. 

Specific star formation rates (SFR per unit mass; 
Guzman et al. 1997) have been obtained using these SFR 
values and the stellar masses given by the highest probabil- 
ity solution cluster in Table [j. The mean SFR and specific 
SFR for SB-like, HII-like, and whole sample are given in 
Table |^ (see also Figure ^d) . 

The specific SFR vs. stellar mass diagram is shown in 
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Figure |^ (see Guzman et al. 1997) . In panel-a we show the 
stellar masses and star formation rates per unit mass for 
three reference samples. We have included the sample of 
Kennicutt (1983, K83 hereafter), taking the Ha and B- 
band luminosities given by K83 a nd th e stellar mass-to- 
light ratios of Faber & Gallagher (1979). In addition, the 
sample of HH-galaxies of Telles (1995) was included, after 



converting virial masses to stellar masses using a correction 
of 0.6 dex (Gallego et al. 1999, in preparation) and assum- 
i ng th e H/3-to-Ha luminosity ratios used by Guzman et al. 
( 1997 ) fo r this sample. Masses and specific SFRs for the 
Calzetti ( 1997b ) sample are also shown. In this case, stel- 
lar masses were inferred subtracting the HI mass from the 
dynam ical mass measured. The SFR values for the Calzetti 
(|l997b|) sample were obtained from their Br7 luminosities 
assuming L HQ = 102.8xL B r 7 (Osterbrook 1989, for T e =10 4 K 
and n e =100cm -3 ). Finally, the dwarf irregular galaxy GR8 



(Reaves 1956) was included. Its Ha luminosity was ob- 
tained from t he H /3 luminosity given by Gallagher, Hunter 
& Bushouse (1989), assuming Lu a / Xh/3=2.86, and its stel- 



l ar ma ss, 3.2x10 Mq, from Carignan, Beaulieu & Freeman 



(1990) . The limits in the Ha luminosity function of Gallego 
et al. fll995j ) 10 40 4 -10 42 - 8 ergs" 1 , are also drawn. 

Figure M shows that the UCM sample clearly repre- 
sents a bridge between relaxed spir al gal axies and the most 
extreme HII galaxies from Telles (1995), that is, Sp— >SB- 
like— >HII-fefce— > HII g alaxies. In fact, some of the HII galax- 
ies from Telles (1995) have very similar properties to those 
of the less massive HH-Zife UCM galaxies, mainly DHIIH 
and B CD s pectroscopic types, very rare in our sample (see 
section 5.1). In addition, most of the SBN type UCM galax- 



ies seem to be normal late-type spirals with enhanced star 
formation. This star formation enhancement is about a fac- 
tor of three, and is due to the ongoing nuclear starburst. 
Thus, the range in specific SFR spanned by the population 
of the star-forming galaxies that dominate the SFR in the 
local universe is (10-10 3 ) x 10" 11 yr" 1 , from the local relaxed 
spirals to the most extreme HII galaxies. In fact, this range 
i s not very different from that obtained by Guzman et al. 
(1997) for a sample of intermediate/high-z compact galaxies 
from the HDF. The high s pecific SFR region, where the HII 
galaxies from Telles ( 1995 ) are placed, is not very well cov- 
ered by our sample due to the scarcity of very low-luminosity 
objects, basically DHIIH and BCD spectroscopic type galax- 
ies, relative to the UCM whole sample. 



dictions was carried out using a maximum-likelihood esti- 
mator in combination with Monte Carlo simulations which 
take into account the observational uncertainties. Our main 
results are: 

(i) The star-forming galaxies in the UCM sample (used to 
determine the SFR density of the local universe), show typ- 
ical burst strengths of about 2 per cent and stellar masses 
of 5xl0 10 Mq. The current star formation in these galaxies 
is taking place in discrete star formation events rather than 
in a continuous fashion. If this is typical of the past star- 
formation history in the galaxies, many of such star forma- 
tion events would be necessary to build up their stellar mass. 
However, our observations provide very little information on 
star-formation episodes that took place before the current 
one. 

(ii) We have identified two separate classes of star- 
forming galaxies in the UCM sample: SB-like and BH-like 
galaxies. Within the HII- like class the DHIIH and BCD spec- 
troscopic type galaxies, i.e. dwarfs, constitute the most ex- 
treme case. The mean burst strength deduced for the SB-like 
galaxies is about a 25 per cent lower than for the dwarf Bll- 
like galaxies. The average stellar mass is an order of mag- 
nitude larger in the former than in the latter. The SB-like 
galaxies are relatively massive galaxies where the current 
star formation episode is a minor event in the build up of 
their stellar masses, while BU-like galaxies are less massive 
systems in which the present star formation could dominate 
in some cases their observed properties and contributes to a 
greater extent to their stellar population. 

(iii) Because of the low burst strengths inferred, A'-band 
luminosity is dominated by the old stellar populations, and 
the Tf-band stellar mass-to-light ratio is almost the same 
(within ~20 per cent) for all the galaxies. Thus, the if-band 
luminosity is a very good estimator of the stellar mass for 
typical star-forming galaxies. 

(iv) The average SFR of the galaxies is log(SFR)~1.5, 
with the SFR expressed in Mq yr -1 , and it is similar for the 
SB-like and the BII-Hke galaxies. However, since the latter 
are typically less massive, their specific SFR (SFR per unit 
stellar mass) is significantly larger, in a factor 2.3, than that 
of the former. 

(v) The UCM galaxies represent a bridge in specific SFR 
between relaxed spirals and extreme HII galaxies. The range 
in specific star formation rate spanned by the local star- 



forming galaxies, (10-10 )x 10 



yr 



" 1 , is very similar to 



that observed in higher redshift objects. 



6 SUMMARY 

Using new nIR observations and published optical data for 
67 galaxies from the Universidad Complutense de Madrid 
(UCM) survey, we have derived the main properties of their 
star-forming events and underlying stellar populations. This 
sample represents about 35 per cent of the UCM galaxies 
covering the whole range of absolute magnitudes, Ha lumi- 
nosities and equivalent widths spanned by the survey. Burst 
strengths and ages, stellar masses, stellar mass-to-light ra- 
tios and, to a certain extent, metallicities, have been ob- 
tained by comparing the observed r — J and J — K colours, 
K-band magnitudes, and Ha equivalent widths and lumi- 
nosities with those predicted by evolutionary synthesis mod- 
els. The comparison of the observations with the model pre- 



APPENDIX: ANALYSIS OF THE SPACE OF 
SOLUTIONS 

In this appendix we briefly describe the analysis performed 
onto the (tlog b, log Z) space of solutions. For each galaxy 
we have 10 3 (£,log &,log Z) points that correspond to the 10 3 
points generated in the (r — J, J — if,2.5xlog EW(Ha)) space 
using a Monte Carlo simulation method. 

The mean values (< t >,<log 6>,<log Z>) are primary 
indicators of the best (t,log fe,log Z) solution, and their stan- 
dard deviations {o t ,(Ti os b,<7\og z) could be taken as estima- 
tors of the deviation of the data. However, due to the well 
known age-metallicity and age-burst strength degeneracies 
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Figure 8. Hypothetical distribution of (t, log b, log Z) solutions. 
PCA1, PCA2, PCA3 are the principal components for each of 
the three solution clusters. 



these standard deviations are not representative of the dis- 
tribution of these solutions in the (t,log fe,log Z) space. For- 
tunately, these degenaracies do not span the whole range 
in age, burst strength and metallicity given by the mod- 
els, being relatively well constrained by the (r — J, J — 
AT,2.5xlogEW(Ha)) data. 

Therefore, we have studied the clustering of the 
(t,log fo,log Z) solutions for each individual galaxy. We have 
used for this analysis a single link age hierarchical clustering 
method (Murtagh & Heck 1987). First, (1) we determine 



the distances between every couple of solutions, which rep- 
resents a total of N x (AT — l)/2 dissimilarities (=distances), 
being N the number of solutions. The dissimilarity between 
the elements j and k, dj t h, is defined as 



d: 



E 



(5) 



The matrix of dissimilarities is known as dendogram. Then, 
(2) we find the smallest dissimilarity, ck k . These points, i 
and k, (3) are therefore agglomerated and replaced with a 
new point, i U k, and the dissimilarities updated such that, 
for all objects j ^ i, k, 



d 



iUk,j 



{dij,dk,j} 



(6) 



Then, (4) the dissimilarities ck,j and d kl j, for all j, are 
deleted, as these are no longer used. Finally, we return to 
step (1) after reducing the dimension of the dissimilarities 
matrix and the number of clusters. Finally, we recover the 
last three clusters of solutions. 

The clustering pattern obtained is basically produced 
by the discretization in metallicity of the original BC96 evo- 
lutionary synthesis models. Now, we analyze the solutions 
within each solution cluster. In this case, the discretizations 
in burst strength and age are comparable and a Principal 
Component Analysis (PCA hereafter) is the most suitable 
choice (0.04 dex in burst strength and ~0.05dex in age). 

The PCA basically determines, in a R n data array, the 
set of n orthogonal axes that better reproduces our data 
distribution. The first new axis, i.e. the principal component, 



will try to go as close as possible through all the data points, 
describing the larger fraction of the data variance. Figure ^ 
shows the principal component for each of the three solution 
clusters of a hypothetical (i,log 6, log Z) distribution. 

Formally, following the PCA (see, e.g., Morrison 1976), 
(1) we construct the variance-covariance and the correla- 
tion matrix of the sample, being the (j, k) th term of these 
matrixes, respectively, 



1 71 

- y2( r v -rj)(r ik -r k ) 



n ^ 



(nj - rj)(r ik - r k ) 



Sj 



where 



(7) 
(8) 

(9) 



Then, (2) solving the eigenvector equation, pu=Xu, we ob- 
tain the eigenvalues and eigenvectors of the correlation ma- 
trix. The ratio between an eigenvalue and the sum of all 
the eigenvalues, Xj/},.^ Xj, gives us the contribution of 
the new axis, determined by the corresponding eigenvector, 
to the total data variance. Therefore, the eigenvector with 
higher eigenvalue is the principal component and will indi- 
cate which is the dominant degeneracy inside each solution 
cluster. 



ACKNOWLEDGMENTS 

This work is based on observations obtained at the Lick Ob- 
servatory, operated by the University of California and on 
observations collected at the German- Spanish Astronomical 
Centre, Calar Alto, Spain, operated by the Max-Planck In- 
stitute fur Astronomie (MPIE), Heidelberg, jointly with the 
Spanish Commission for Astronomy. It is also partly based 
on observations made with the Isaac Newton Telescope op- 
erated on the island of La Palma by the Royal Greenwich 
Observatory in the Spanish Observatorio del Roque de los 
Muchachos of the Instituto de AstroHsica de Canarias. 

A. Gil de Paz thanks the Institute of Astronomy of the 
University of Cambridge for all the facilities and support 
during his stay there. J. Gallego and A. Gil de Paz ac- 
knowledge the invitation, hospitality and facilities provided 
during the 3rd Guillermo Haro Workshop, included in the 
Guillermo Haro Programme at the INAOE (Mexico). We 
also thank C.E. Garcfa-Dabo, C. Sanchez Contreras and R. 
Guzman for stimulating conversations and the referee Dr. 
M. Edmunds for his useful comments and suggestions. A. 
Gil de Paz acknowledges the receipt of a Formation del Pro- 
fesorado Universitario fellowship from the Spanish Ministry 
of Education. A. Aragon-Salamanca acknowledges generous 
financial support from the Royal Society. A. Alonso-Herrero 
was supported by NASA on grant NAG 5-3042. This re- 
search was also supported by the Spanish Programa Secto- 
rial de Promocion General del Conocimiento under grants 
PB96-0610 and PB96-0645. This work was partially carried 
out under the auspices of EAR A, a European Association for 
Research in Astronomy, and the TMR Network on Galaxy 



© 1999 RAS, MNRAS 000, 0-|l| 



18 A. Gil de Paz et al. 

Formation and Evolution funded by the European Commis- 
sion. 



REFERENCES 

Abraham R. G., Ellis R. S., Fabian A. C, Tanvir N. R., Glaze- 
brook K., 1999, MNRAS, 303, 641 

Alonso O., Zamorano J., Rego M., Gallego J., 1995, A&AS, 113, 
399 

Alonso O., Garci'a-Dabo C. E., Zamorano J., Gallego J., Rego M., 

1999, ApJS, 122, 415 
Alonso-Herrero A., Aragon-Salamanca A., Zamorano J., Rego M., 

1996, MNRAS, 278, 417 (AH96) 
Aragon-Salamanca A., Ellis R. S., Couch W. J., Carter D., 1993, 

MNRAS, 262, 764 
Armand C, Baluteau J.-P., Joubert M., Gry C, Cox P., 1996, 

A&A, 306, 593 
Balzano V. A., 1983, ApJ, 268, 602 
Barger A. J. et al., 1998, Nature, 394, 24 

Baugh C. M., Cole S., Frenk C. S., Lacey C. G., 1998, ApJ, 498, 
504 

Bernasconi P. A., Maeder A., 1996, A&A, 307, 829 

Bershady M. A., Haynes M. P., Giovanelli R., Andersen D. R., 

1999, in Galaxy Dynamics, ASP conference series, 146, 235 
Bland-Hawthorn J., Maloney P. R., 1997, PASA, 14, 59 
Brocklehurst M., 1971, MNRAS, 153, 471 
Bruzual A. G., Chariot S., 1993, ApJ, 405, 538 (BC93) 
Calzetti D., Kinney A. L., Storchi-Bergmann, T., 1996, ApJ, 458, 

132 

Calzetti D., 1997a, in Waller H. et al., eds, AIP Conf. Proceedings, 
Vol.408, The Ultraviolet Universe at Low and High Redshift: 
Probing the Progress of Galaxy Evolution. American Institute 
of Physics, New York, p. 403 

Calzetti D., 1997b, AJ, 113, 162 

Cardelli J. A., Clayton G. C, Mathis J. S., 1989, ApJ, 345, 245 
Carignan C, Beaulieu S., Freeman K. C, 1990, AJ, 99, 178 
Chariot S., 1998, in Benvenuti P. et al., eds., The Next Genera- 
tion Space Telescope: Science Drivers and Technological Chal- 
lenges, 34th Liege Astrophysics Colloquium, ESA-SP, 429, 135 
Connolly A. J., Szalay A. S., Dickinson M., Subbarao M. U, 

Brunner R. J., 1997, ApJ, 476, Lll 
Courteau S., A compilation of the UKIRT Faint JHK Standars, 

NOAO, 1995 
Coziol R., 1996, A&A, 309, 345 

Driver S. P., Windhorst R. A., Griffiths R. E., 1995, ApJ, 453, 48 
Elias J. H., Frogel J. A., Matthews K., Neugebauer G., 1982, AJ, 
87, 1029 

Ellis R. S., 1997, ARA&A, 35, 389 
Faber S. M., Gallagher J. S., 1979, ARA&A, 17, 135 
Fernie J. D., 1983, PASP, 95, 782 
Fioc M., Rocca-Volmerange B., 1999, A&A, 344, 393 
Gallagher J. S., Hunter D. A., Bushouse H., 1989, AJ, 97, 700 
Gallego J., 1995, Ph.D. thesis, Universidad Complutense de 
Madrid 

Gallego J., Zamorano J., Aragon-Salamanca A., Rego M., 1995, 
ApJ, 455, LI 

Gallego J., Zamorano J., Rego M., Alonso O., Vitores A. G., 1996, 

A&AS, 120, 323 (GAL96) 
Gallego J., Zamorano J., Rego M., Vitores A. G., 1997, ApJ, 475, 

502 

Gallego J., 1998, in Thuan T. X., Balkowski C, Cayette V., Tran 
Tranh Van J., eds., Dwarf Galaxies and Cosmology. Proceed- 
ings of the XVIIIth Moriond astrophysics meeting, Editions 
Frontieres. Gif-sur-Yvette, France 

Glazebrook K., Blake C, Economou F., Lilly S., Colless M., 1999, 
MNRAS, 306, 843 



Guzman R., Koo D. C, Faber S. M., Illingworth G. D., Takamiya 
M., Kron R. G., Bershady M. A., 1996, ApJ, 460, L5 

Guzman R. et al., 1997, ApJ, 489, 559 

Hammer F. et al., 1997, ApJ, 481, 49 

Heraudeau P., Simicn F., 1997, A&A, 326, 897 

Herbst T. M., Beckwith S. V., Birk C, Hippler S., McCaughrean 
M. J., Mannucci F., Wolf J., 1993, in Fowler A. M., ed., Proc. 
SPIE, Vol. 1946, Infrared Detectors and Instrumentation. The 
International Society for Optical Engineering, p. 605 

Hu E. M., Cowie L., McMahon R. G., 1998, ApJ, 502, 99 

Hughes D. et al., 1998, Nature, 394, 241 

Kennicutt R. C, 1983, ApJ, 272, 54 

Kennicutt R. C, 1992, ApJ, 388, 310 

Kent S. M., 1985, PASP, 97, 165 

Kriiger H., Fritze-v. Alvensleben U., Loose H.-H., 1995, A&A, 
303, 41 

Leitherer C, Ferguson H. C, Heckman T. M., Lowenthal J. D., 

1995, ApJ, 454, L19 
Lilly S. J., Le Fevre O., Crampton D., Hammer F., Tresse L., 

1995, ApJ, 455, 60 
Lilly S. J. et al., 1998, ApJ, 500, 75 

Lowenthal J. D., Koo D. C, Guzman R., Gallego J., Phillips A. 
C, Faber S. M., Vogt N. P., Illingworth CD., Gronwall C, 

1997, ApJ, 481, 673 

Madau P., 1997, in Holt S. S. & Mundy L. G., eds., AIP Conf. 

Proceedings, Vol.393, Star Formation Near and Far. American 

Institute of Physics, New York, p. 481 
Madau P., Ferguson H. C, Dickinson M. E., Giavalisco M., Steidel 

C. C, Fruchter A., 1996, MNRAS, 283, 1388 
Madau P., Dickinson M. E., Pozzeti L., 1998, ApJ, 498, 106 
Mathis J. S., 1990, ARA&A, 28, 37 

Misch A., Gilmore K., Rank D., 1995, UCO-Lick Technical Re- 
ports, 77, Sta. Cruz, California 

Mobasher B., Sharpies R. M., Ellis R. S., 1993, MNRAS, 263, 560 

Morrison D. F., 1976, Multivariate Statistical Methods, McGraw- 
Hill Book Co., Singapore 

Murtagh F., Heck A., 1987, Multivariate Data Analysis, D. Reidel 
Publishing Co., Dordrecht, Holland 

Osterbrock D. E., 1989, Astrophysics of Gaseous Nebulae and 
Active Galactic Nuclei. University Science Books, Mill Valley, 
California 

Pei Y. C, Fall S. M., 1995, ApJ, 454, 69 
Peletier R. F., Balcells M., 1996, AJ, 111, 2238 
Poggianti B. M., 1997, A&AS, 122, 399 
Reaves C, 1956, AJ, 61, 69 

Rowan-Robinson M. et al., 1997, MNRAS, 289, 490 
Salzer J. J., MacAlpine G. M., Boroson T. A., 1989, ApJS, 70, 
479 

Scalo J. M., 1986, Fund. Cosmic Phys., 11, 1 

Somerville R. S., Primack J. R., Faber S. M., 1999, MNRAS, in 
press 

Steidel C. C, Giavaliso M., Pettini M., Dickinson M., Adelberger 
K. L., 1996, ApJ, 462, 17 

Steidel C. C, Adelberger K. L., Giavaliso M., Dickinson M., Pet- 
tini M., 1999, ApJ, 519, 1 

Storchi-Bergmann T., Calzetti D., Kinney A. L., 1994, ApJ, 429, 
572 

Telles E., 1995, Ph.D. thesis, Cambridge University 
Thuan T. X., Gunn J. E., 1976, PASP, 88, 543 
Tresse L., Maddox S. J., 1998, ApJ, 495, 691 

Treyer M. A., Ellis R. S., Milliard B., Donas J., Bridges T. J., 

1998, MNRAS, 300, 303 

Vitores A. G., 1994, Ph.D. thesis, Universidad Complutense de 
Madrid 

Vitores A. G., Zamorano J., Rego M., Alonso O., Gallego J., 

1996a, A&AS, 118, 7 
Vitores A. G., Zamorano J., Rego M., Gallego J., Alonso O., 

1996b, A&AS, 120, 385 



© 1999 RAS, MNRAS 000, 



Star formation properties of local galaxies 

Wainscoat R. J., Cowie L. L., 1992, AJ, 103, 332 
Williams R. E. et al., 1996, AJ, 112, 1335 
Worthey G., 1994, ApJS, 95, 107 

Zamorano J., Rego M., Gallego J., Vitores A. G., Gonzalcz- 

Riestra R., Rodriguez G., 1994, ApJS, 95, 387 
Zamorano J., Gallego J., Rego M., Vitores A. G., Alonso O., 1996, 

ApJS, 462, 1 

This paper has been produced using the Royal Astronomical 
Society/Blackwell Science M^X style file. 



© 1999 RAS, MNRAS 000, 0-|l| 



